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INTRODUCTION 


by  CDK  Keith  N<  Sargent)  USN 


From  the  standpoint  of  National  Defense  the  most  vexing  problem 
facing  the  Surveillance  and  Detection  Planners  is  the  combination  of 
physical  limitations  imposed  by  line  of  sight  propagation  from  a  tan¬ 
gent  point  on  a  sphere  and  the  requirement  to  detect  in  the  total 
volume  surrounding  the  sphere^  With  the  advent  of  higher  velocity 
vehicles  together  with  the  resultant  shortening  of  available  reaction 
time  and  the  exponential  increase  in  aero-space  volume  to  be  searched 
as  orbital  dimensions  grow,  the  need  for  a  breakthrough  in  surveil¬ 
lance  capability  became  more  and  more  urgent. 

*  As  a  result,  when  Dr»  S.  Fred  Singer  and  his  University  of 
Maryland  group  approached  the  Bureau  of  Ships  with  his  hypothesis 
regarding  plasma  sheath  anomalies  in  space,  the  Detection  and 
Countermeasures  Branch  of  the  Warfare  Systems  Research  and  Develop¬ 
ment  Division  was  receptive  to  his  thinking.  If  his  hypothesis 
proved  out,  it  was  recognized  that  target  enhancement  could  be 
achieved  that  would  permit  tracking  of  satellites  at  ranges  where 
skin-tracking  was  impracticable,  if  not  indeed  impossible.  It  was 
further  recognized  that  plasma  sheath  detachments  due  to  coulomb 
drag  could  provide  track  signatures  to  detectors  before  the  vehicle 
itself  rose  above  the  radar  horizon. 

Therefore,  a  modest  effort  was  financed  through  the  Electro- Jet  \ 
Corporation  to  permit  Dr,  Singer  and  his  group  to  pursue  his  hypothesisX 
This  report  is  the  result  of  their  efforts.  While  experimental  proof  ' 
has  not  yet  been  obtained,  an  experiment  has  been  designed  and  will 
be  performed  to  this  end. 

Because  of  the  promise  for  a  unique  satellite  and  missile  early 
warning  capability  inrougu  me  i^roua-iauger  erru-.t ,  j  t  is  censidared 
that  this  technical  report  should  be  snared  with  tne  scientific  military 
community  in  advance  of  the  experimental,  proof.  It  is  anticipatad 
that  a  subsequent  technicaj.  report  will  be  issued  upon  completion  and 
evaluation  of  the  experiment  which  is  planned  for  the  near  tuLura, 


tExperltaental  Studies  of  the  Kraus  Effect 
by  RIchari  T.  Bettinger 

The  following  is  based  on  my  trip  to  the  Ohio  State  Observatory  on  June 
21,  1961 j  Or.  Kraus'  paper  (Ref.  i)  and  Dr.  Villard's  paper  (Ref.  2). 

Critical  Resume  of  Observations 

The  Ohio  Stata  woall  iMs  (acluded  both  radar  and  C.W.  type  observations, 

the  former  extending  from  early  19^  through  late  I96O.  The  radar  equipment 

utilized  for  the  most  part  consisted  of  a  2^  KW  peak  po««er  transmltt^  oper* 

ating  at  a  pulse  repetition  rata  of  30  pulses  per  sec.,  a  pulse  width  of  one 

millisecond,  and  a  drequaacy  range  from  17~I9  me.  The  antennas  utilized  included 

V- 

everything  from  1/2  wave  dipole  to  a  fixed  20  DB  parabolic  reflector.  The  most 
cennonly  used  receiver  band  width  was  b  kc  with  values  varying  from  1  •>  8  kc. 

The  Creft-Villard  experiments  at  Stanford  used  two  systems.  The  first 
consisted  of  a  taw  radar  transmitter  operating  with  1  millisecond  pulses  at 

a  repetition  rate  of  10  pulses  per  sec.  Frequencies  of  from  IO-3O  me  were  utilized  and 
aioderately  high  gain  antennas  (fixed  rhoaiblcsi  were  in  generel  aa^leyed.  The 
Stanford  group  also  included  a  phase  coherence  aaaaureamntf  l.a.  phase  ef  the 
returning  signal  was  eampared  to  that  of  the  generating  system  and  the  phase 
difference  recorded.  No  statement  concerning  receiver  band  width  is  madei 
however,  it  is  presumed  to  be  of  the  order  of  2  kc. 

Except  for  the  additional  phase  measurement  at  Stanford,  there  seems  to 
be  no  fundamental  difference  in  the  experimental  techniques  aaployed  by  the  two 
groups.  The  Stanford  group  had  available  a  considerably  wider  frequency  range, 
however,  they  state  that  the  system  was  adjusted  for  maxiiaun  return  dependent 
t^n  the  particular  ionospheric  conditions  that  exist  at  the  time.  This  would 
ieply  that  fraquenctes  of  the  order  of  1?  me  dominated,  and  the  average  frequency 
used  should  net  differ  significantly  from  of  thia  State. 
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The  data  taking  techniques  and  interpretation  of  the  data  of  the  two  groups 
differed  significantly  however.  The  Stanford  group  considered  two  satellites  and 
set  up  to  make  measurements  cniy  on  those  occasions  when  these  two  satellites 
would  pass  through  the  beams  of  one  of  their  available  antennas.  The  Ohio  State 
groupj  however^,  tended  to  set  up  their  equipment  and  allow  it  to  run  continuously^ 
and  then  observe  the  data  obtained  for  “interesting  or  unusual  events".  When  one 
of  these  was  detected  they  attempted  to  the  best  of  their  ability  to  relate  it  to 
a  satellite  pass.  This  is  a  difficult  process  because  of  the  relatively  large 
number  of  satellites  in  orbit  and  the  dearth  of  readily  available  accurate  data 
on  all  of  these.  This  process  of  locating  all  satellites  at  the  time  of  this 
event  was  obviously  time  consuming  and  has  prevented  a  large  number  of  events 
from  being  examined.  The  Stanford  group,  on  the  other  hand,  defined  their  interval 
of  Interest  as  the  period  in  which  the  satellite  was  between  the  half  power  points 
and  counted  all  anomalies  in  that  time  as  possible  events.  The  events  were  classi¬ 
fied  to  the  best  of  the  reader's  ability  in  terms  of  the  probability  that  it  was  a 
satel 1 i te- induced  or  satellite-related  event.  The  details  of  the  rating  system 
will  not  be  gone  into  here,  howeve*-  the  techniques  employed  seem  reasonable 
and  under  the  circumstances  as  good  as  can  be  done  with  this  type  of  approach. 

The  Ohio  State  group  employed  still  another  technique  for  the  measurement 
of  these  reflection  anomalies  wnich  very  simply  consisted  of  a  receiving  antenna 
and  receiver  set  to  WWV  at  10,  or  20  me.  The  antennas  again  varied  considerably 
but  In  general  were  relatively  low  gain.  The  most  important  aspect  of  this  experi¬ 
ment  was  that  the  receiver  had  relatively  poor  long  term  stability  so  that  the 
center  frequency  over  a  period  of  time  tended  to  drift  up  to  ^  kc  from  its  prior 
setting.  Receiver  band  widths  were  of  the  order  of  2  or  ^  kc.  Although  originally 
set  directly  on  WWV,  after  running  overnight  they  were  very  often  plus  or  minus 
5  kc  from  the  original  setting.  In  this  off-tune  position  a  number  of  events 
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which  Kraus  has  described  as  satellite  related,  have  been  detected.  A  high 
velocity  plasma  cloud,,  (or  reflecting  surface  of  some  type)  will  Doppler  shift 
the  signal  (depending  upon  orientation),  and  place  the  reflected  signal  outside 
of  the  receiver  capabilities  of  either  of  the  radar  experiments  previously  described. 
However  because  of  the  accidental  drift  in  the  Kraus  CW  experiment  they  would 
presuaubly  measure  these  Doppler  shifted  signals  with  a  considerable  enhancement 
of  signal  to  no!.,.-  ^  :i. .  Widening  the  receiver  band  width  sufficiently  to  include 

a  widely  shifted  Doppler  signal  would  have  degraded  the  signal  to  noise  ratio 
badly,  possibly  making  the  signals  undetectable.  The  evaluation  techniques  again 
are  similar.  An  event  is  noted  and  an  attempt  is  made  to  correlate  it  to  a 
satellite  passage,  if  some  reasonable  correlation  exists  then  other  factors  are 
introduced  as  possible  variables.  A  large  number  of  these  events  have  been 
reported  in  detail  by  Kraus. 

Discussion 

The  following  observations  seem  germane: 

1.  Little  or  no  attention  has  been  ascribed  to  the  ability  or  inability  of  the 
antenna  systems  used  to  localise  the  source  of  the  reflection?  we  are  simply 
trading  off  reflecting  area  versus  antenna  orientation,  i.e.  the  same  response 
is  obtained  from  a  s!de  lobe  which  ss  10  db  down  and  looking  at  a  tenfold  larger 
reflecting  surface.  In  genera'  the  techniques  employed  to  date  provide  no 
localization  In  space  beyond  a  very  rough  indication. 

2.  When  considering  the  reflection  caused  by  the  local  ionization  in  the 
vicinity  of  a  satellite,  one  in  general  would  expect  little  or  no  results. 

If  the  satellite  is  charged  negatively  then  the  charge  sheath  consists  of  a 
patitiva  foe  sheath  whose  thickness  is  in  general  only  a  faw  dabye  lengths 
and  whase  projected  area  will  scarcely  exceed  that  of  tho  physical  area  af 
satellites.  In  any  case,  radio  propagation  is  relatively  insensitive  to  the 


ion  charge  sheath,  which  merely  represents  a  small  hoi j  in  Llio  electron  distri¬ 
bution  which  viil!  tend  to  s.taitcr  bet  not  reflect.  The  immediate  wake  of  the 
satellite  does  represent  an  electron  enhanced  area,  more  precisely  an  ion 
evaluated  area,  but  one  would  not  expect  considerably  enhanced  electron  den¬ 
sities.  In  short,  one  would  not,  except  under  exceptional  circumstances, 
expect  any  excessive  ionizaiion  to  be  present  in  the  immediate  vicinity  of 
the  satellite.  There  will  be  a  trail  of  ionization  turbulence  in  the  wake 
of  the  satellites  which  I  would  off  hand  assume  to  be  relatively  small.  In 
no  case  has  anyone  reported,  that  I  know  of,  direct  reflections  from  ionization 
trails  related  to  satellite  passages.  In  general,  Kraus'  work  which  has  been 
the  motivation  of  these  general  studies  has  indicated  reflections  from  anomalies 
which  tend  to  originate  at  or  in  the  general  vicinity  of  a  satellite  orbit,  but 
are  in  general  moving  in  directions  not  associated  with  the  orbit  and  in  general 
at  very  high  velocities.  We  would  conclude,  therefore,  that  an  experiment 
(such  as  the  Stanford  one)  which  strictly  localizes  the  immediate  vicinity  of 
the  satellites  would  probably  not  be  fruitful. 

3.  The  events  reported  by  Kraus'  are  sufficiently  few  that  the  effects,  if  in 
fact  satellite  related,  are  relatively  rare  and  presumably  related  to  fairly 
specialized  conditions.  Hence,  any  statistical  study  must  have  a  relatively 
large  source  in  order  to  establish  the  presence  of  correlation. 

At  this  stage  of  the  game  alt  the  work  that  has  gone  into  this  area 
still  seems  highly  inconclusive.  Kraus'  measurements  may  in  fact  indicate 
satellite  correlated  local  areas  of  ionization,  however,  this  is  certainly 
not  conclusive  based  on  his  evidence.  This  work  is  in  no  way  prejudiced  by 
the  Stanford  report  since  they  allowed  themselves  almost  no  possibilities  of 
detecting  the  type  of  phenomenon  reported  by  Kraus.  If  in  fact  these  high 
velocity  areas  of  enhanced  ionization  do  exist,  then  the  relatively  narrow 
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band  widths  of  the  receivers  used  in  the  radar  studies  almost  preclude  the 
possibility  of  detection.  (The  exception  of  course  is  a  cloud  moving  with 
zero  radial  velocity.) 

The  Stanford  work,  on  the  other  hand,  seems  to  establish  that  there  are 
no  directly  associated  ion  clouds  with  reflective  areas  larger  than  I  or  2 
orders  of  laagnltude  greater  than  the  satellite  physical  area  directly  assoc* 
lated  with  the  satellite  while  moving  in  orbit.  The  numbers  are  still 
sufficiently  small  so  that  the  statistical  conclusion  is  subject  to  reason¬ 
ably  large  errors,  however  no  large  scale  systematic  effect  Is  present.  The 
result  is  of  course  not  surprising. 

The  nature  of  the  phenomnon  with  which  we  are  dealing  is  best  described 
by  Kraus;  *^st  of  our  observations  Indicate  a  transient  or  sporadic  phenomenon 
in  which  detectable  ionization  appears  in  the  general  vicinity  of  the  satellite, 
but  not  necessarily  very  close  to  it.  9n  fact,  fast  moving  echoes,  ^Ith  velocities 
of  from  30  to  130  kilonMtert  per  sec  which  appear  to  originate  In  the  setelllte 
seem  to  be  frequently  involved”.  Kraus"  data  seems  to  Indicate  that  these 
clouds  do  not  propagate  parallel  or  perpendicular  to  the  magnetic  field,  but 
In  general  at  some  skew  angle.  We  might  summarize  the  germane  characteristics 
of  such  a  phenomenon  In  terms  of  the  experimental  parameters  of  prime  Interest. 

1.  Effective  area  as  a  function  of  frequency. 

2.  Space-time  localization  with  the  highest  possible  accuracy. 

3*  Effective  cloud  velocity. 

4.  Trail  reflectivity  and  persistence,  if  it  exists. 

The  experimental  problems  associated  with  obtaining  definitive  Information 
on  these  parameters  are  obviously  trasmndous,  however  if  a  serious  attampt  to 
describe  this  type  of  phenomenon  is  to  be  made,  these  parameters  must  be 
attacked  simultaneously.  The  velocity  of  the  clouds  is  a  very  important 
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polnt;  a  cloud  moving  with  a  radial  velocity  of  100  kilomatars  par  aae  will 
hava  a  Dappler  shift  of  12  kc,  hence  it  wiil  in  general  be  outside  tha  response 
band  of  the  radars  currently  in  use,  hence  undetectable.  The  motion  relative 
to  the  earth's  magnetic  field  is  presumably  not  random  so  that  the  motion 
relative  to  the  observing  station  will  probably  show  preferential  relative 
angles  and  a  limited  range  of  radial  velocities. 
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PLASMA  SMEATH  ANO  SCREENING  AROUND  A  RAPIDLY  MOVING  BODY 


E.  H.  Walker 

ABSTRACT 

The  distribution  of  the  potential  and  the  charge  density  is  derived 
quite  generally  for  a  stationary  charged  sphere  and  for  a  charged  body 
moving  rapidly  through  a  plasma.  Previous  treatments  were  restricted  to  the 
cases  where  either  the  body's  potential  was  small,  being  limited  to  less  than 
five  times  the  plasma  temperature  (i,e.  I8  ;6r  5kT/e) ,  or  the  body  was  small 
compared  to  the  Debye  length  (  i  ,e.  a  <  (kT/i+tm^e*)  ),  or  the  body  was  moving 
slowly  compared  to  the  ion  velocity. 


PLASMA  SHEATH  AM'>  SrRF.Ef^lKG  AROUND  A  RAPIDLY  MOVING  BODY 


1.  INTRODUCTION 

Several  workers^^*' have  investigated  the  nature  of  screening 
of  bodies  of  various  geometries  that  are  embedded  in  a  plasma.  This  question 
is  particularly  important  in  considering  such  problems  as  the  significance  of 
electron  and  ion  density  measurements  In  the  laboratory  and  in  space,  the 
charging  of  probes  and  satell  ites,  and  the  production  of  a  wake  by  such 
bodies . 

The  limitations  that  beset  the  existing  calculations  ares 

(a)  The  assumption  of  lota!  thermodynamic  equilibrium  is  very  often 
unsatisfactory.  In  many  practical  screening  problems,  under  laboratory  or 
space  conditions,  the  mean  free  path  of  particles  is  found  to  be  large  com¬ 
pared  with  the  Debye  screering  length. 

(b)  The  use  of  a  barometric  type  of  formula  for  the  density  of 
screening  particles  is  generally  not  valid.  Although  this  is  equivalent  to 
the  previous  assumption,  it  is  sometimes  used  even  though  calculations  are 
performed  that  take  account  of  the  motions  of  the  particles.  This  is  done 
either  because  calculations  are  limited  to  small  bodies  or  because  of  a 
failure  to  distinguish  between  ronact ret i ng  and  accreting  orbits, 

( c)  The  assumption  that  the  ions  do  not.  respond  quickly  enough  to  the 
potential  field  and  thus  maintain  their  amoient  density  is  a  simplifying 
assumption  that  is  not  sat isfactory.  (I)  in  the  case  of  the  screening  of 

a  stationary  body,  there  is  no  distinction  to  be  made  between  the  screening 
of  a  positively  charged  body  and  a  negatively  charged  body.  (2)  It  must  be 
noted  that  as  e  particle  with  energy  kT  moves  into  a  region  where  the  potential 
is  greater  than  kT,  It  must  be  subjected  to  considerable  acceleration.  It  will 
be  found  that  for  satellite  velocities,  considerable  deflections  do  occur. 
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( d)  The  calculations  generally  apply  only  to  either  small  potentials  or 
to  small  bodies  However,  since  potentials  on  satellites  probably  reach  values 
of  50  kT,  and  since  the  bodies'  dimensions  are  usually  many  times  the  Debye 
screening  distance,  these  limitations  are  unsatisfactory  for  most  problems. 

We  wish  to  obtain  solutions  of  the  screening  of  a  body  which  will  not  be 
subject  to  these  limitations.  In  particular,  we  will  obtain  exact  solutions 
to  the  pfLo'em  of  screening  for  the  case  of  a  stationary  sphere  embedded  in  a 
plasma  and  also  we  will  obtain  the  screening  for  a  rapidly  moving  body. 

2.  SELF-CONSISTENT  CALCULATION  OF  THE  SCREENING  OF  A  STATIONARY  CHARGED  SPHERE. 

We  wish  to  obtain  solutions  for  the  problem  of  the  screening  of  a  stationary 
charged  sphere  that  is  embedded  in  a  plasma  with  a  mean  free  path  for  its  con¬ 
stituents  that  is  long  compared  to  the  Debye  length  h.  The  general  equations 
that  will  be  used  here  have  been  derived  by  Opik^.  A  brief  resume  of  the 
essentia!  details  is  presented  below. 

In  the.  neighborhood  of  any  point  in  a  spherically  symmetric  potential 
field  we  may  express  the  potential  in  the  form 

0  =  Kr’*^  (1) 

where  K  and  n  must  be  evaluated  in  that  neighborhood.  It  will  be  convenient 
for  the  calculations  to  use  this  type  of  expression. 

The  conservation  of  angula*"  momentum  and  energy  may  be  written  in  terms 
of  the  particle  velocity,  v,  at  r,  the  impact  parameter,  q,  and  the  velocity 
of  the  particle  at  infinity,  u  (See  Fig,  1). 


rv  sinct  ■  qu 

(2) 

^/amv®  ■  0  +  ^/amu® 

(5) 

where  a  is  the  angle  between  ^  and  £,  and  0  is  the  potential  energy  of  the 
particle  relative  to  infinity  and  is  positive  for  an  attractive  force  and 


negative  for  a  repulsive  force.  If  we  write 
U  “ 

we  have  on  combining  the  angu'ar  momentum  and  energy  equations 

sin^a  >  (4) 

This  expression  can  be  used  to  distinguish  between  two  types  of  particle 
motions.  The  first  type  is  called  periastron  motion.  For  this  type  as  r 
decreases  from  infinity,  there  always  exists  a  value  of  r  for  which 

q*  -  r^  (l+ei/U)  (5) 

so  that 

s  i  a  “  ! 

Under  this  condition  the  space  inside  th's  value  of  r  is  not  accessible  to 
the  particle.  For  a  repulsive  ootentiai  this  condition  can  always  be 
satisfied. 

If  we  now  substitute  Eq,  (!)  into  Eq,  (5)  we  have 

sin*a  =“  - - - q^/F  (6) 


When  n<2,  the  denor, i natcr,  F,  always  app’-oaches  zero  as  r  goes  to  zero. 

Uowe-i/er,  it  ”>”2  t^'C'c  is  a  point  for  which  sir.^Ct  reaches  a  maximum 
value.  This  value  of  i  s  ob’-ained  by  setting  the  derivati-e  of  F  equal  to 

zero 

-  0 


dF 


giving 

2  +  20/U  +  r21'/U  «  C 
which  gives  for  the  potential  at 


2U 


c  (n  -2' 
'  c 


(8) 


(9) 
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The  value  of  r  is 


2U  J 


(10) 


where  the  subscripts  refer  to  the  values  of  the  quantities  at  corresponding 
value  of  r^.  Thus,  inside  r^  the  particle  spirals  toward  the  origin.  Here 
the  orbit  is  described  as  pericritical .  For  a  given  value  of  q,  the  maximum 
value  of  a  is  obtained  from  Eq.  (4) 


( s  i  n®a) 


max 


(„  -2)  2/" 


(11) 


We  note  that  n^  and  K^.  depend  on  U  but  on  nothing  else.  Hence  (sin®  )^^^ 
varies  directly  as  q®.  Evidently,  the  limiting  condition  for  entry  into  the 
periastron  region  is 

(IS) 

Thus  Eq.  (11)  sets  a  lower  limit  on  q^  for  the  target  radius  of  periastron 
type  orbits  for  the  case  n>  2 

'IS  * 


n 

c 

- 


1/n. 


(13) 


(nc-2)^"c 


or 


q  =  r  /  n  /(n  -2'. 
^p  c  /  c  c  • 


(14) 


Particles  satisfying  Eq,  (13)  will  enter  a  periastron  and  emerge  from 
it  in  a  pseudo-hyperbolic  trajectory.  Those  with  q q^,  however,  will  always 
make  an  angle  of  less  than  9^°  with  the  radius  and  spiral  inwards  indefinitely. 
These  orbits  are  of  the  pericritical  type. 

We  are  now  in  a  position  to  calculate  the  accretion  of  particles  through 
a  spherical  surface.  The  number  of  particles  passing  through  a  target  ring  q 
to  q  +  dq  is 
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dA^  =  Strq  dq  uNq 

where  Nq  is  the  ambient  density.  Now  from  Eq,  (4)  we  have 

q  dq  =  (  '  +  l5/U)  sin  Q(  cos  Qt  dcK  (16) 


Thus 


dAo(  =•  27rr^uNo(l  +  iS/U)  sinOt  costK  d  Ct  (17) 

The  inteyral  of  this  for  the  periastron  case  for  which  all  values  of  Cl  are 
allowed  is 

Ap^  =  2Trr^uNo(  I  +  0/U)  j  sinCl  cosCl  dd 

0 

=  irr^u  No( !  +  0/U)  ( 1 8) 

Now  fo  ■  pericritica!  interaction,  when  (S  ^  0,  n>  2  and  r  C  r^,  the  accretion 

is  constant,  involving  all  the  particles  that  have  reached  the  pericritlcal 

surface,  r  .  Thus,  the  accretion  has  the  same  value  for  all  r  <  r 
c  c 

Ap^  -  Trr^®  uN^(  I  +  0^/U)  (19) 

This  result  is  also  obtained  by  using  the  limits  0  to  Cl  in  the  integration 
of  Eq.  (18)  . 

The  particle  densitv  can  now  be  obtained  in  terms  of  the  potential 
energy  from  the  above  accretion  formula.  Thus,  we  have  the  density  contributed 
by  particles  moving  at  angle  between  U  and  U  +  dCl 


dN  =  dA/4Tri^v  cosCl  -  Va  No(  l+lfUl  sinCl  dot  (20) 

For  periastron  orbits  an  integration  of  this  yields 

N  »  V2QNo(1  +  I5/U)  sinCt  da=  V2(iNo(l  +  l9/U)  (21) 

P®  J 

0 

The  factor  4  Is  the  return  factor.  If  a  particle  that  passes  through  the 
spherical  surface  at  r  is  elastically  reflected  by  the  body  so  as  to  return 
and  again  pass  through  the  surface  at  r,  the  particle  orbit  must  be  counted 
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twice  and,  thus,  the  return  factor  will  have  a  value  2  for  this  particle.  If 
the  particle  is  absorbed  and  does  not  return,  the  value  is  1  for  this  orbit. 

The  return  factor  in  Eq.  (21)  is  an  average  value  of  the  return  factor  taken 
over  all  particles.  If  the  collision  is  inelastic,  a  more  involved  calculation 
may  become  necessary.  It  will  be  found,  however,  that  inelastic  collisions 
are  not  important  since  the  reflection  coefficient  for  both  ions  and  electrons 
is  smal 1®< Furthermore,  although  secondary  electrons  are  produced, 
they  ordinarily  can  be  neglected;  a  calculation  that  includes  their  effect 
will  be  made  in  Section  4. 


where  N_  and  for  a  negatively  charged  sphere  are  respectively,  the  density 
of  electrons  and  the  density  of  the  ions  times  Z.,  the  number  of  charges  per 
ion.  For  a  positively  charged  sphere,  is  the  density  of  electrons  while 
N_  is  the  density  of  the  ions  times  Z. . 

Equation  (23)  can  be  written  in  non-dimensional  form  by  defining 

if  "^/kT;  p  “  r/h  (24) 

where  k  is  the  Boltzmann  constant,  T  is  the  temperature  of  the  plasma  and  h 
is  the  Oebye  screening  distance  as  given  by 


-T‘ 


where  e  is  the  electronic  charge  and  Nq  is  the  ambient  density.  Thus  Eq.  (2'ji) 
becomes 


N.  N+  0^1 
1,11  - - Z.  ~  =1. 


(26) 


■  No  No  p 

and  for  a  negatively  charged  body  embedded  in  a  plasma  of  singly  ch-i-i'  c!  ions 
the  densities  are 


where 


w 


N_/No  -  V2a(!  - 

N^/No  =  VaQi  1  ■*-  w 


1, 


1/n 

when  n  <  2  or  (ln-21^/3)  <  I 


Va 

in  al 1  other  cases 


(27) 

(28) 


(29) 


The  value  of  n  must  be  evaluated  at  each  point  in  the  integration  since 
the  potential  does  not  follow  a  simple  power  law  over  a  wide  range  of  r. 

It  will  be  noticed  that  a  Maxwellian  velocity  distribution  has  not  been 
included  in  the  derivation  of  these  equations.  Expressions  have  been  obtained 
which  take  account  of  the  Maxwe>lian  velocity  distribution^,  but  because  of 
the  added  diffic^-ltY  of  rhe  cal cul at  ions,  these  equations  have  not  yet  been 
empi oyed , 


3  RESULTS  OF  THE  SCREENING  CALCULATIONS 

The  integration  of  Eq.  (26)  using  Eqs.  (27),  (28)  and  (29)  for  densities 
of  the  screening  particles  was  obtained  by  numerical  methods.  The  self-con¬ 
sistent  field  calculation  proceeded  from  a  position  po,  for  which  a  small 
value  of  the  potential  \|ro  was  specified,  toward  the  origin  (rather  than 
starting  from  the  surface  of  the  charged  body)  The  formula  derived  for  the 
screening  under  local  thermodynamic  equilibrium  provided  an  approximate  value 
for  s 
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♦  '  “  -  { 1  +  p)*/p  (50) 

any  error  in  the  initial  conditions  will  be  automatically  damped  out  with  a 
proper  choice  for  the  integration  interval  ^  p.  This  procedure  produced  curves 
that  are  independent  of  the  radius  of  the  charged  sphere  provided  either  that 
few  electrons  reach  the  surface  of  the  body  or  that  the  value  of  Q  is  2.  We 
have  set  the  return  factor  for  the  electrons  and  for  the  ions  to  be  equal  in 
this  calculation.  It  can  be  included  as  an  additional  parameter  in  f  by 
setting 

=  (51) 

It  is  this  function  which  is  plotted  in  the  curves. 

In  Fig.  2  we  have  plotted  log  against  log  p,  to  obtain  a  family  of 
curves.  A  curve  connecting  the  pericritical  surfaces  is  also  shown.  It  will 
be  found  from  these  curves  that  the  potential  falls  off  much  more  slowly  than 
in  the  case  of  local  thermodynamic  equilibrium.  In  Fig.  Ja,  b  we  plot  the  ratio 
h*/h  against  p,  where  h*  is  the  value  of  the  Debye  length 

necessary  in  the  Debye  type  screening  formula  0^  to  make 
0p(h*,r)  =  0  (h,r)  and  0p  (h*, r)  =  0  ;h^r) 

in  our  calculation  for  a  given  point-  A  set  of  curves  are  given 

which  correspond  to  the  set  of  curves  in  Fig.  2.  Thus  the  Debye  type 
formulation  gives  a  much  stronger  screening  than  that  presented  here. 

The  most  interesting  curves,  however,  are  those  for 

M_j./Mo,  N_/Mo  plotted  against  p,  as  shown  in  Figs.  4a,  b,  c,  d,  e,  f,  g.  In 
these  curves,  the  density  of  the  electrons  (assuming  the  body  to  be  negatively 
charged)  falls  rapidly  to  zero,  whereas  the  density  of  the  ions  first  rises 
slightly,  then  falls  rapidly  to  quite  low  values  and  finally  rises  rapidly  as  the 
screening  particles  converge  toward  the  origin. 
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4,  THE  bCKEENING  OF  A  CHARGED  BODY  WHEN  IS  VARIABLE 

In  experiments  performed  by  Fowler  and  Farn.sworth^^  it  has  been  found 
that  the  reflection  coefficient  for  low  energy  electrons  striking  a  metallic 
surface  iS;  for  our  purposes,  small;  its  value  is  of  the  order  of  5°/o' 
Hagstrum®' has  found  that  ions  striking  a  metal  surface  will  be  neutral¬ 
ized  by  Auger  charge  exchange  with  an  efficiency  of  nearly  100°/o«  A 
would  therefore  seem  that  it  is  not  satisfactory  to  set  (j,  equal  to  2  for 
particles  which  strike  the  surface  of  the  body.  We  will  derive  an  c -.ij  ress  i  on 
for  the  return  factor  for  the  ions,  Q,^,  and  we  will  assume  a  sufficiently 
large  negative  potential  for  the  body  to  allow  us  to  use  a  value  of  2  for 
the  electron  return  factor,  The  value  of  can  be  obtained  from  the 

expression  for  the  accretion.  If  A^  represents  the  accretion  through  a 
surface  of  radius  p  and  represents  the  accretion  through  the  pericritical 
surface,  the  value  of  the  ionic  return  factor  outside  of  the  pericritical 
surface  is 


Ac  o 

A 

P 


PC^(  I  +  ®/3  'I'c) 

P®TT+ 


(>2) 


where  the  subscripts  refer  to  the  values  of  the  quantities  on  the  pericritical 
surface.. 

The  use  of  two  different  expressions  for  the  ionic  and  electronic  return 
factors  modifies  the  screening  most  strongly  for  small  values  of  If  we 
assume  ♦  to  be  small,  the  Poisson  equation  can  be  written  as 

-  VaQ,  (1  -  -  Vai^  (1  + 


If  we  take  p  to  be  large  enough  then  on  the  pericritical  surface  as  obtained 
from  our  previous  cal  cut ar ions,  will  be  less  than  1.  Thus  the  first  derivative 
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term  in  Eq.  (33)  small  (these  conditions  are  met  if  p  is  about  10  or 

more)  and  we  can  neglect  it.  We  will  then  have  the  condition 


a  -  (L  a  +  (L  , ,, 


(J*) 


which  must  be  satisfied  in  order  for  the  screening  equation  to  be  valid. 
Substituting  for  on  the  left  side  of  the  inequality  we  obtain 


U\  >  ..  t  I  T  -/a  ivc) 

( 1  +  “/3 i+i ) 


(I  %l»c\) 


(55) 


If  we  further  reduce  this  using 
I 

2 

we  obtain^  approximately: 

lil  >  ®/a(f)^ 

We  also  have  the  condition  that  if  no  screening  at  all  occurs 

this  shows  that  'Ir  must  vary  as 

t=Kr*^i;  l<:.n<2 


(36) 


(37) 

(38) 

(39) 


for  small  values  of  Furthermore,  this  result  will  hold  in  a  modified 

form  for  the  case  of  local  thermodynamic  equilibrium. 

The  per i cr i t i cal  surface  as  indicated  by  Eqs.  (37)  and  (38)  must  be 
given  by 

Pj.^/*/3  kl'  P  ('^O) 

If  0^  ■  2,  then  Eq.  (3^)  is  always  satisfied,  and  for  the  region  of  small 
if  an  Integral  of  the  equation  can  be  obtained  of  the  form 

k;)"^  e"^  (41) 

It  might  appear  that  if  one  chooses  p  to  be  large  enough,  then  the  difference 
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between  and  <1_  could  be  neglected.  By  the  substitution  of  Eq.  (32)  into 
the  Poisson  equation  we  can  stiow  what  prevents  this.  We  have 


^1 

P 


^  - 


(42) 


for  small  \|f.  Thus,  if  we  choose 

‘^/a  1  ilr  1  (45) 

we  would  obtain  a  potential  of  the  form  given  in  Eq.  (45).  However,  p^,  p, 
and  \|f  cannot  be  chosen  independently.  Furthermore,  if 

p^2/2p® -jt  2/3  j^l  (44) 

then  since  +  varies  as  p  the  condition  in  Eq.  (43)  cannot  be  satisfied  for 
any  other  choice  of  ’I'  and  p.  It  appears  that  the  potential  does  not  vary 
according  to  Eq.  (4l)  over  any  portion  of  the  curve  for  which  is  small  and 
0^  is  given  by  Eq.  (32) . 


Experiments  performed  by  Hagstrum®^ have  shown  that  as  the  result  of 
the  neutralization  of  ions  striking  a  metal  surface,  secondary  electrons  will 
be  produced  with  an  efficiency  of  about  25°/o‘  These  electrons  have  not  been 
included  in  the  above  considerations.  Their  effect  may  be  easily  included 
into  the  calculations  by  introducing  into  the  Poisson  equation  the  term  N^/Nq 
where 


,  iTPr^  (  I  •» 

4p^  /  ‘^/3(  to  ■*  If)  '' 


(45) 


where  7  is  the  efficiency  of  secondary  electron  p'-oduction,  and  to 
potential  at  the  surface  of  the  charged  body  (or  alternatively  at  the  surface 
for  which  the  secondary  electron  energy  is  zero). 

In  Fig.  2,  we  have  plotted  a  curve  to  show  the  effect  of  this  term  on  the 
screening.  Nevertheless,  the  effects  of  the  secondary  electron  term  upon  the 
previous  results  will  be  small  if  the  surface  of  the  body  has  a  high  potential. 
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3.  SELF-CONSISTENT  FIELD  CALCULATION  FOK  A  MOVING  CHARGED  BODY. 

The  calculation  of  the  screening  of  a  moving  charged  body,  except  for 
cases  of  extremely  low  velocities,  is  much  more  complicated  than  the  calculations 
for  a  stationary  body.  The  accretion  and  density  formulae  obtained  in  the  second 
section  depend  heavily  upon  the  spherical  symmetry  of  the  problem.  This  symmetry 
allows  us  to  integrate  over  a  and  q  without  consideration  of  variations  of  the 
integrand  over  an  equ ipotent ial  surface  and  without  consideration  of  individual 
particle  orbi ts. 

For  the  calculation  of  the  screening  of  a  moving  body,  it  will  be  necessary 
to  follow  individual  particle  orbits  in  order  to  calculate  density  contributions. 
These  density  contributions  will  then  be  used  to  calculate  the  protential  in  the 
neighborhood  of  the  particle  position,  and  then  the  potential  is  used  to  con¬ 
tinue  the  orbit  calculations,  etc. 

In  order  to  avoid  the  necessity  ofiterations  in  the  calculations  we 
assume  to  be  unity  for  all  particles.  Thus  a  particle  will  be  artificially 
removed  from  the  calculation  when  its  orbit  becomes  parallel  to  the  local 
equipotential  surface.  It  should  be  noted  from  the  solution  of  the  screening 
of  a  stationary  sphere  that  if  mo>  th®  initial  position  for  the  calculation 
is  chosen  to  be  large,  then  the  pericritical  surface  occurs  almost  immediately. 
Thus,  very  few  particles  will  be  excluded  in  an  artificial  way.  The  actual 
computations  will  proceed  as  follows. 

First,  the  position  of  a  new  equipotential  surface,  0^,  is  found,  making 

use  of  the  orbit  density  as  obtained  for  the  0^  position  of  the  particle  from 

a  calculation  of  orbits  from  the  0  ,  to  the  0  surface.  The  calculation  of 

n-i  n 

the  potential  follows  lines  which  we  will  call  0  lines,  perpendicular  to  the 
equipotential  surfaces  (actually,  these  are  the  electric  field  lines).  Cal¬ 
culations  are  made  from  the  0^_j  to  the  0^  surface  by  making  use  of  the  appropriate 
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value  of  the  line  element  perpendicular  to  the  t  surface.  For  this  cal¬ 
culation,  the  position,  the  direction  of  the  surface  mormai  and  the  value  of  0' 

n 

must  be  obtained. 

Secondly,  the  orbit  density  is  found.  Since  the  new  potential  surface 
has  been  obtained,  the  orbits  can  be  calculated  starting  from  the  last  0  surface 
and  proceeding  to  the  new  0  surface.  In  order  to  obtain  densities  from  these 
orbits,  however,  it  is  necessary  to  also  obtain  the  position  of  another  particle 
orbit  that  Is  nearby.  We  will  obtain  later  an  expression  for  the  separation 
of  the  nearby  particle.  A  weight  factor  Is  carried  In  the  calculations  for  each 
particle  orbit  which  represents  the  spacing  of  the  particles  at  infinity.  This 
weight  factor  will  be  used  in  making  averages  over  several  orbit  densities. 

Thus,  if  we  initially  start  with  one  region  on  the  0  surface  (actually  the 
quantity  we  are  considering  is  the  initial  Impact  pareMters  rather  than  the 
initial  angle)  thickly  sown  with  test  particles,  whereas,  another  region  is 
sparbeiy  sown,  the  densities  contributed  by  the  particles  in  the  sparse  region 
must  count  more  than  the  other  particles.  Each  orbit,  along  with  its  companion, 
contains  enough  information  to  give  a  value  of  the  density  at  the  position  of 
the  particle,  not  Just  a  density  contribution.  However,  it  is  necessary  to 
average  over  several  densities  when  several  orbits  pass  near  a  0  line  position. 

6.  DENSITY  OF  COLO  SCREENING  PARTICLES  IN  AN  AXIALLY  SYMMETRIC  POTENTIAL  FIELD 
WITH  VELOCin  U. 

Assume  that  we  have  an  axially  symmetric  charged  body  moving  in  a  cold 
plasma,  in  the  center  of  mass,  we  can  view  the  body  as  stationary  while  ions 
move  parallel  to  the  axis  of  symmetry  with  a  velocity  u  and  impact  parameters 
raprasantad  by  q.  The  accretion  of  particles  between  q  and  q  ^  (sea  Fig  .  3) 

Is 


M  -  Sirq^  N«u 
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The  density  N  can  be  related  to  the  flux  J  and  the  velocity  of  the 
particles  at  some  point  v  by 

N  -  J/v  (47) 

Now  the  flux  is  given  In  terms  of  the  component  of  the  axial  surface  element 
perpendicular  to  v,A  by 

J  -  M/ASj.  (48) 

Since 

ASj.  *  *1'*  ^  ^  (49) 

where  O  |s  the  angle  between  the  velocity  vector  and  the  radius  vector, 
we  have  tor  the  density 


,  45^  ,  iJ  nr.«5  cosoT  r*;,  ♦  (/u)^»  cos  a  * 

Now,  in  terms  of  the  angle  subtended  by  we  can  write 

8ln«o-^  (51) 

t<here  the  subscripts  refer  to  the  initial  values  of  these  quantities  on  the 
initial  potential  surface  as  specified  by  the  boundary  conditions  that  we 
choose.  The  differential  Is 


cos9o  ^0 


iOa 

ro 


(52) 


so 

H  .  No  ro^  sinQy  cosQq  /Sgo - 

r*(l  +  0/U)  cosCl  sin©  ZiP 

It  will  be  noticed  that  N  goes  to  infinity  as  a  approaches  the  value  90**. 
This  is  due  to  the  neglect  In  the  above  derivation  of  the  curvature  of  the  sur¬ 
face  element  t&  which  Is  subtended  by  the  angle  We  correct  this  by  obtain¬ 
ing  tStx.  by  nn  integration. 

By  axial  sysiMtry  we  can  write  for  ibSa.(aee  Fig.  6  ) . 

ASj.  •  2irr  sin©  (54) 

where  AJl  is  the  projection  of  the  width  of  the  surface  ring  element  perpendicular 
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to  the  particle  orbit,  if  is  divided  Into  elements  d^Pwe  have  for  these 

dAj-  rd^j  sin  (90®  “  a  +  Z!0,)  (55) 

integration  yields 

Aj}"v|^cosas!n^+sina(l  -  cos  j  ( 56) 

Therefore,  we  can  write 

ASj_  ■  2irr*  sin®  |^cos  C(sin^‘*‘sina(l  ■=  cos  £0)  j  (57) 

rather  than 


2Tr*  sin®  cosa  £0 
We  have  for  N 


(58) 


No  ro*  sin®o  cos  ®o  ^o 

m - i.  — -  -  -  --  — [ — 1  1 - -I - -  1 

r^l  +  0/U)  sin®  [  cos  a  sIn/K®  'f'  since  (I  cos  ^)]  (59) 

A  further  correction  must  be  made  for  the  sin®  factor.  As  0  approaches 
zero  N  approaches  infinity.  This  effect,  which  is  quite  accurate  under  the 
assumption  of  a  cold  plasma,  is,  nevertheless,  unsatisfactory  and  arises 
from  our  neglect  of  the  therstal  velocities  of  the  screening  ions,  initially 
the  test  particles  in  velocity  space  fill  a  sphere  of  radius  v.^.  where 


located  at  a  distance  u  frotri  the  origin.  All  particles  with  velocity  vectors 
initially  lying  in  the  plane  of  the  calculations  will  have  for  their  minimum 
angular  approach  to  the  axis,  9  .  ,  a  value  zero.  All  test  particles  with  v- 
lying  out  of  this  plane  will  pass  with  some  larger  value  of  If  i  is  the 

total  intograted  path  of  the  particles,  then  the  average  of  Is 

obtained  with  the  aid  of  fig,  7s 

1'  ■  I  ?  <«'» 

If  we  defina  "K  and  X'  to  be  the  distances  and  j['  in  terms  of  labye  lengths 
X  -  ]/»•  {eg) 


(6>) 


Then 
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ir.h  \.;1iich 


Pq*  sin©o  cosSo  ^0 

p*(l  +  U/U)  cos(l©  -  900|  -  ^  ) 

cosa  s  1  n^s  I  nCi^  1  -cos/i0) j 

-!■-  ;l:.!  ;.!!  •. 

iniJuui". 

"In' 

cos  (le  -  so”!  -  5„,„) 

The  expression  for  the  density  then  becomes 

N  -  No 


7.  CALCULATION  OF  FROM  £Bo 

»■ 

In  our  calculation  of  particle  densities  w*  wish  to  follow  one  tesC 
particle  orbit  and  obtain  by  appropriate  computations  the  value  of  the 
density  at  the  points  along  this  orbit.  Thus  far  we  have  obtained  an 
expression  for  the  density  involving  the  value  of  C(,  d,  do>  Po>  P> 

£00  and  £0.  |f  we  know  the  particle  orbUty.  All  of  these  quantities  may  b^ 
Obtaineil  except  for  the  quantities  £0o  and  £0.  These  presuppose  the  existence 
of  two  orbits  and  these  quantities  designate  the  angular  separation  of  the 
particle  positions  when  they  are  on  equi potential  surfaces.  The  first  of 
these,  £00,  designates  the  initial  separation  of  the  particles. 

We  now  wish  to  obtain  £0  by  assuming  the  two  orbits  to  be  extremely 
close.  Thus,  we  will  be  able  to  find  an  expression  giving  the  change  in 
£0  as  we  proceed  along  only  one  orbit. 

An  express ion^for  d©  in  terms  of  dr,  J^,  the  angular  momentum,  U,  the 
initial  kenetic  energy,  v,  the  distance  to  the  particle,  and0(r},  the 
potential  energy  in  a  central  field.  Locally  we  can  assume  our  potential 
function  to  be  a  central  field.  Thus,  we  will  write 
_ dy _ 


(64) 


d©  - 


mv»  /|(U  +  «(r) 


(6^) 


ffluq 


Wt) 


(6$) 


where 
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Now  ^9  the  change  In  the  angular  separation  between  the  two  particle 
positions  as  shown  In  Fig ^  8  does  not  depend  explicitly  upon  r,  0,  nor  dr, 
since  these  quantities  on  the  same  equi potential  surface  will  be  the  same 
for  both  particles.  Also,  since  initially  the  two  particles  had  the  same 
energy,  this  Is  not  the  variable  that  determines  the  variation  in  d9.  The 
angular  momentum,  however,  is  not  the  same  for  the  two  particles.  Initially  It 
Is  given  by 

0  *  01  u<l  *  xMir  slnOo  (6B) 

and  £p  la 


-  mu  r  sin  (Oq  ♦  •  £o  *  Aiio 


Now 


a(< 


r  |^sin(ao  *  ^o)  "  alocto  j  ■  muq  cotOo 


(69) 

(70) 


or  more  generally  between  two  surfaces 

m  >!  r  cosa  (Tl) 

It  should  be  noted  that  the  assumption  of  the  conservation  of  momentum  is  not 
necessary  to  this  derivation  since  the  particles  are  close  together.  Variations 
in  the  angular  momentum  will  only  contribute  higher  order  terms  to 


Now,  let  us  differentiate  d0  with  respect  toi.  This  gives 

m  sin®  a  1 


Add  ■  d9 


!  + 


’U+0 


(72) 


and  upon  substitution  for  aJ! 

Ad9  -  dO  cotO  [l  ♦  2  tan*d  j  (73) 

Where  dO  is  the  angular  awvement  of  the  particles  In  moving  from  one  potential 
surface  to  the  next  and  A9  Is  the  separation  of  the  particles. 

These  formulae  along  with  those  used  to  follow  individual  particles  from 
surface  to  surface  and  to  follow  the  0  lines  give  the  equipment  for  a  self-con¬ 
sistent  field  calculation  of  both  the  screened  electric  field  and  the  density 
of  screening  ions. 


-18- 


8.  RESULTS  OF  THE  SCREENING  CALCULATION  FOR  A  TYPICAL  PROBLEM 

The  results  of  one  set  of  calculations  are  shown  in  Fig.  9  and  11.  The 

calculations  are  made  for  a  body  moving  at  a  velocity  u  ■  7  km/sec  through  a 

plasma  of  0*  ions  and  electrons,  corresponding  to  the  altitude  range  of  500 

to  1500  km  in  the  earth's  atmosphere.  The  temperature  of  the  plasma  is  taken 

as  T  “  1500°C.  The  radial  distance  from  the  center  of  the  body  to  the  edge 

of  the  screening  sheath  is  set  at  I80  Debye  lengths.  Of  course,  the  shape  and 

size  of  the  body  itself  will  be  taken  to  be  that  of  the  appropriate  equipotential 

surface  (to  be  determined  from  the  solution  of  this  problem).  The  values  of 

are  given  in  terms  of  the  kinetic  energy  of  the  ambient  ions  as  measured  in 

% 

the  coordinate  system  moving  with  the  satellite,  i.e.  m.u®  and  not  in 
terms  of  kT;  the  relation  between  the  two  is; 

♦(0/kT)  -  .0318  . 

The  electrons  have  been  included  in  this  calculation  by  merely  setting 

N_  -  Noe"®^*"^  (74) 

where  N_  is  the  density  at  a  potential  t  in  terms  of  the  ambient  density  Nq- 
This  expression  is  satisfactory  for  our  purpose  and  can  be  shown  to  be  exact 
under  all  conditions  involving  an  exclusively  repulsive  field,  even  where  the 
magnetic  field  is  important. 

It  will  be  noticed  that  Fig.  9  includes  results  only  on  the  front  side 
of  the  body.  This  is  due  to  the  fact  that,  as  yet,  satisfactory  initial 
conditions  have  not  been  calculated  for  the  back  side  of  the  body.  The  cal¬ 
culations  given  actually  involved  following  51  electric  field  lines  and  75 
particle  pai rs. 

In  Fig.  9  we  give  a  plot  of  the  equipotential  surfaces,  the  equidenslty 
surfaces  of  the  ions  and  six  simple  trajectories  of  particles  moving  In  the 
field.  These  values  are  displayed  using  a  plot  with  the  radius  as  the  ordinate 
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( in  units  of  Debye  lengths)  and  the  angle  0  (in  radians)  plotted  along  the 
abscissa  (the  corresponding  distance  is  also  provided).  Figure  10  indicates 
the  mapping  involved  in  this  type  of  plot.  It  should  be  noticed  that  the 
scale  of  the  ordinate  is  smaller  by  a  factor  of  100  as  compared  to  the 
abscissa.  Thus  the  trajectories  of  the  ions  shown  in  Fig.  9  are  actually 
quite  curved.  On  the  other  hand,  the  equ i potent ial  and  equidensity  surfaces 
are  only  slightly  aspherical. 

In  Fig.  11  we  show  plots  of  the  density  of  ions  and  potential  as  a 
function  of  the  radial  distance  for  three  different  angles  to  the  direction 
of  motion. 
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Diagram  showing  the  parameters  involved  in  the  motion  of  a  screening 
particle. 

if  plotted  again  p^,  with  pericritical  surface  indicated  and  one  curve 
with  secondary  electrons.  The  curves  a,  b,  c,  d,  e,  f,  and  g  represent 
the  family  of  curves  that  give  the  potential  as  a  function  of  the 
radius  vector.  For  example,  assume  a  body  with  radius  of  23  Debye 
lengths  has  a  potential  of  15  kT  volts.  This  would  be  point  A  and 
thus  the  potential  would  fall  off  according  to  curve  e.  The  peri¬ 
critical  surface  would  occur  at  a  distance  of  3^  h  and  potential  at 
0.105  kT,  If  we  include  the  effects  of  secondary  electrons  and  charge 
this  sphere  to  a  potential  20  kT  then  this  will  put  us  at  point  c  so 
that  the  curve  e'  represents  the  variation  of  ilr  vs.  p. 

A  plot  of  hVh,  the  ratio  of  the  effective  Debye  length  to  the  actual 
value,  vs.  the  nondimens ional  distance  p.  These  curves  correspond 
to  those  of  Fig.  2  as  indicated  by  the  letter  labels. 

Plot  of  N^/No  and  h'  /No  against  p  for  the  potential  curve  a. 

Plot  of  N^No  and  N  /No  against  p  for  the  potential  curve  b. 

Plot  of  N^/Nq  and  N  /Nq  against  p  for  the  potential  curve  c. 

Plot  of  N^/No  and  N  /Nq  against  p  for  the  potential  curve  d. 

Plot  of  N^/No  .^nd  N  /No  against  p  for  the  potential  curve  e. 

Plot  of  N_^/No  and  N  /Nq  against  p  for  the  potential  curve  f. 

Plot  of  N^/No  and  N  /Nq  against  p  for  the  potential  curve  g. 

The  initial  parameters  involved  in  the  self-consistent  calculation  of 
screening  for  a  moving  sphere. 

Illustration  of  the  parameters  used  to  calculate  the  correction  to 
A  Sj^  for  smal  1  a. 


Fig.  7. 


Fig.  6. 


Fig.  9. 


Fig.  10. 


Fig.  II. 


Illustration  of  the  parameters  used  to  calculate  0^.^  for  the  correction 
of  sin0. 

Illustration  of  the  change  Ad©  in  the  angular  separation  of  the 
companion  particle. 

A  plot  of  the  equ ipotential  surfaces,  equidensity  surfaces  and  six 
sample  particle  trajectories  for  the  front  side  of  a  satellite.  The 
body  is  assumed  to  be  moving  at  7  km/sec  through  a  plasma  of  0^  ions 
and  electrons.  The  temperature  of  the  plasma  is  taken  to  be  1500°C. 

The  ordinate  is  the  radial  distance  in  Debye  lengths  and  the  abcissa 
gives  the  angle  in  radians. 

Figure  to  show  the  mapping  employed  in  Fig.  8.  A  portion  of  a  polar 
plot  has  been  mapped  into  a  cartesian  coordinate  system. 

The  variation  in  the  nondimens ional  potential  (&/kT  and  the  normalized 
density  n  along  the  radius  vector  for  several  representative  angles: 
0=0,  1.005,  1.445  radians. 


Fig.  3  • 


Fig.  3  b 


Figure 
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PLASMA  COMPRESSION  EFFECTS  PRODUCED  BY  SPACE  VEHICLES 
IN  A  MAGNETO  IONIC  MEDIUM 

S.F.  Singer  and  E.  H.  Walker 


ABSTRACT 

\  novel  phanomenon  is  discussed  which  leads  to  the  creation  of  large 
clouds  of  compressed  plasma^.  It  is  based  on  a  theoretically  predicted 
mechanism  which  is  caused  by  the  interaction  of  an  electrically  charged 
body  with  the  surrounding  ambient  plasma,  in  the  presence  of  a  magnetic 
field.  The  cloud  consists  of  ions  which  have  been  scattered  into  quasi-trapped 
orbits.  Their  injection  and  removal  are  studied  in  detail.  Under  non-steady 
State  conditions  the  cloud  may  detach  itself  from  the  body  and  give  rise  to 
"ghQSts".  It  Is  quite  possible  that  this  phenomenon  is  the  cause  of  the 
unusua!  "bursts"  and  radar  echoes  observed  to  be  associated  with  passes  of 
earth  satel I i tes . 


I  .  Statement  of  the  Problem 


If  a  body  carrying  an  electric  charge  moves  rapidly  through  an  ionized 
gas  pervaded  by  a  magnetic  field,  ions  will  be  deflected  by  the  electrostatic 
field  of  the  body  and  will  remain  quasi-trapped  by  the  magnetic  field  until 
removed  by  scattering.  We  wish  to  consider  the  extent  to  which  a  charged  body 
can  produce  a  buildup  of  the  ion  density  along  or  near  its  path. 

if  we  consider  the  ion  velocities  to  be  much  less  than  the  velocity  of 
the  charged  body,  i.e.  if  the  plasma  is  cool,  then  in  a  coordinate  system 
moving  with  the  charged  body  the  incident  ions  (field  particles)  all  move  with 
a  velocity  v  until  they  are  deflected  through  an  angle  d  (which  depends  on  the 
impact  parameter  p) ,  and  then  either  move  ahead  of  or  behind  the  charged  body 
(Fig.  I).  Because  of  thepresence  of  the  magnetic  field  (assumed  for  the 
moment  to  be  parallel  to  v)  those  ions  (or  test  particles)  that  are  deflected 
through  an  angle  close  to  90®  will  continue  to  move  with  the  deflecting  body 
and  will  remain  in  the  neighborhood  of  the  body  until  scattered. 

The  buildup  of  the  density  distribution  of  the  test  particles  is  countered 
by  four  possible  removal  mechanisms: 

Collision  with  ambient  ions  (field  particles), 
c.  Collision  with  neutral  atoms, 

^  Co'lisicn  with  other  test  particles, 

4.  Collision  with  the  charged  body. 

The  collisions  with  the  neutral  atoms,  because  of  the  low  cross  section, 
will  not  be  as  important  as  collisions  with  the  ambient  ions  at  altitudes  above 
^00  km.  Collisions  with  other  test  particles  will  only  be  important  when  the 
compression  of  the  test  particles  becomes  large;  furthermore,  their  scattering 
will  not  be  qualitatively  different  from  that  of  the  field  particles,  and  there¬ 
fore  they  can  be  Included  by  an  appropriate  choice  of  the  mean  free  path  \  of 
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the  test  particles. 

If  the  potential  of  the  body  is  positive  and  greater  than  kT/Ze,  collisions 
with  the  charged  body  will  not  be  important.  However,  special  attention  must 
be  given  to  the  case  of  a  negative  sphere.  This  will  be  considered  along  with 
the  effects  of  screening  in  Section  9*  Por  an  initial  approach,  we  will  con¬ 
sider  the  case  of  a  low  density  plasma  so  that  we  may  neglect  the  effects  of 
screening. 

2.  Calculations  of  Important  Parameters  for  the  Deflection  of  Ions  by  Charged 
Bod i es 

It  is  important  to  consider  the  orders  of  magnitude  of  several  lengths 
which  enter  into  our  problem.  The  following  equations  are  in  HKS  units. 

We  have  the  following  relation  between  the  impact  parameter,  p,  and  the 
scattering  angle  6 

p  -  f  ■  Po  cot  I  «  9-57  X  lO"^  cot  I  (I) 

where  E  is  the  energy  of  the  ions,  q  the  charge  on  the  body,  V  is  the  potential 
on  a  body  of  radius  a,  Ze  is  the  charge  of  the  ion,  A  is  the  atomic  weight  of 
the  ion  and  po  is  the  90°  impact  parameter. 

The  radius  of  gyration,  p  ,  of  a  particle  in  a  magnetic  field  is  given  by 

p  ■  sin  ©  ■  po  s  i  n®  ■  1 . 045  xl0®^sin©  (2) 


where  m  is  the  mass  of  the  particle,  B  is  the  magnetic  induction  and  Po  is  the 
radius  of  gyration  for  a  test  particle  moving  at  right  angles  to  the  field. 


The  Debye  screening  length  h  is 


T 


where  k  Is  Boltzmann's  constant,  T  is  the  temperature  of  the  plasma,  and  n  is 
the  ambient  electron  density. 
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Table  1  contains  numerical  values  of  these  important  quantities. 

The  mean  free  path  for  energy  loss  and  for  angular  scattering  for  an  ion 
moving  somewhat  faster  than  the  thermal  ions  but  more  slowly  than  the  thermal 
electrons  is  important  in  determining  our  approach  to  the  problem  of  the  build 
up  of  an  ion  cloud  near  a  satellite. 

For  an  ion  which  is  moving  at  a  velocity  of  7  km/sec  through  a  plasma  of 
ionized  oxygen  at  a  temperature  of  1500°K  (for  which  the  root  mean  square 
thermal  velocity  for  the  ions  is  1.53  km/sec  and  for  the  electrons  261  km/sec) 
the  coefficient  of  dynamical  friction  produced  by  the  electrons  divided 

by  <j2w,,>.  for  the  ions  is,  according  to  Spitzer* 


<  ^  .  (■"iV  .055)  .  ^  . 

<  >j  2z.“  m.  G(w/w.)  2G(6.0) 


(^) 


and  for  hydrogen  ions  the  values  is  0.051.  For  this  reason  we  need  not  be 
concerned  with  energy  loss  to  the  electrons.  The  ratios  of  the  deflection 
time  tjj  to  the  energy  exchange  time  t^  and  also  to  the  slowing-down  time  t^, 
again  according  to  Spitzer^,  are 


V'e  •  ■  °-°5T 

for  oxygen  and  1.12  for  hydrogen  while 


(5) 


Vn  •  (5:)  "-oo  (6) 

for  oxygen  and  0.7^  for  hydrogen.  This  shows  that  change  in  the  velocity  of 
the  test  particles  will  occur  approximately  as  fast  as  the  change  in  angle. 

In  the  derivation  of  an  expression  for  the  density  and  extension  of  a  buildup 
of  test  particles  into  a  cloud  moving  with  the  charged  body,  we  will  first 
consider  the  angular  deflection  of  the  test  particles,  assuming  the  velocity 
to  be  constant  and  then  consider  the  effect  of  their  slowing  down. 

The  effective  cross  section  a  for  the  deflection  of  a  test  particle 
by  the  field  particles  can  be  obtained  from  Spitzer*: 


(cgs  units) 


(7) 


in  A 


For  V  ■  7  kffl/sec,  T  ■  1500®K,  n  ■  10®  cm  *  the  mean  free  path  ^  for  a  90® 
deflection  of  an  0^  test  particle  is  900  km;  for  the  mean  free  path  is 


^  km. 

31.  Approximate  Density  of  the  Deflected  Ions 

An  approximate  expression  for  the  number  density  or  concentration  of  the 
deflected  test  particles,  n,  is  obtained  by  dividing  (7rp‘vT)n  ,  the  volume 
swept  out  by  the  charged  body  in  a  time  T  times  the  ambient  ion  (or  electron) 
density,  by  rriZp)^^,  the  volume  containing  the  helical  paths  of  the  ions  for  a 
time  interval  T.  Since  we  have  for  X 

/  ■  vT  cos  fl  (8) 


we  get  for  n 


cos® 


(9) 


Let  us  take  as  an  example  a  spherical  satellite  of  radius  I  meter  charged 
to  a  potential  of  10  volts  moving  with  a  velocity  of  7  km/sec  in  a 
magnetic  field  of  ?  x  10~®  webers/n^.  If  we  take  en  Intermediete  value  for 
6  of  60°  then  n/n  has  a  value  of  10.9  for  ions  but  is  completely  negligible 
for  0*^  ions  since  they  would  be  moving  too  rapidly  to  be  affected  appreciably 
by  the  magnetic  field. 

A  more  refined  expression  than  Eq.  (9)  can  be  obtained.  Let  us  assume 
that  In  an  equilibrium  situation  the  rate  of  injection  of  test  particles 
Into  helical  paths  corresponding  to  a  deflection  angle  0  Is  Iq*  If  these 
test  particles  are  removed  by  collisions,  what  will  bo  their  number  at 
any  time  t  after  the  Injection?  We  have 
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where  the  relative  velocity  between  the  test  particles  and  the  field  particles, 


v^,  is  given  by 


V  -  2v  S  in  ;r 
r  2 


Therefore,  we  have 


-Hgovrt 


(11) 


(12) 


The  rate  of  injection  of  test  particles  between  the  angles  9  and  9  +  d9 
can  be  expressed  in  terms  of  £  (9),  the  scattering  cross  section  for  ions 
deflected  by  the  charged  sphere: 


d9  -  J  1(9)  2-rr  sin9  d9 

where  J  is  the  flux -of  incident  field  particles  given  by 

J  ■  n  V 
e 

£  (9)  is  given  by 


cos  9/2 


so  that 


•  Po  2sin9"s'in^ 


9  COS  9/2  -OeOVrt 

^9  ’  "e"  711^972  « 


(13) 


(14) 


(15) 


(16) 


The  time  can  be  expressed  in  terms  of  the  distance  from  the  center  of  the 
charged  sphere,  x,  by  the  substitution  for  t  of 

t  ■  -  x/( V  cos  9)  (IT) 


If  we  substitute  Eq.  (I6)  into  (l4)  and  integrate  over  x  we  obtain 
1^  ^  _  2qnesin  9/2 

I  -  e  cos  9 


9 


TTPf 


.  cos  9/2 


2»  sin^  9/2 


(18) 


We  can  drop  the  exponent  if  we  take  x  to  be  large.  Integrating  over  all 
angles  from  9o  to  l80^  gives  v  for  the  total  number  of  ions  N 


irpoT 

3  0 


sin 


■a  ®a 


I 


(19) 


If  n^  is  defined  as  the  linear  density  of  ions,  its  average  value 


-6- 


Equations  (2^)  and  (26)  give  us  the  equilibrium  density  of  test  particles 
along  the  path  of  the  charged  body.  We  will  use  these  expressions  to  determine 
the  compression  produced  by  satellites  moving  parallel  to  a  magnetic  field. 

In  Eq.  (6)  we  showed  that  the  mean  time  for  the  slowing  down  of  the  test 
particles  is  almost  the  same  as  the  mean  time  for  the  deflection  of  the  test 
particles  for  both  and  0*.  This  will  mean  that  the  forward  velocity  of  the 
test  particles  will  fall  off  as  an  exponential  function  (or  even  faster) 
about  as  rapidly  as  the  particles  are  removed  by  deflections.  Since  the 
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density  varies  inversely  with  the  velocity  of  the  ions,  this  effect  results 
in  a  sustained  or  in  some  cases  a  growing  density.  The  sweep  back  of  these 
ions  is  also  important.  This  gives  rise  to  a  folding  back  of  the  ion  cloud 
onto  itself.  For  large  clouds  such  that  Po^>  Po  the  sweep  back  ions  will 
contribute  as  greatly  as  the  main  stream  of  ions  but  where  po>>  Po  this 
would  not  be  the  case  except  for  the  fact  that  the  ions  are  slowing  down 
and  reducing  po-  We  should  therefore  add,  {n  either  case^  a  factor  of  2  to 
Eqs.  (25)  and  (26)  to  account  for  the  sweep  back  effects. 

4.  The  Diffusion  Equation 

Let  f(x,9)  dx  dll  represent  the  number  of  ions  in  a  plane  perpendicular 
to  X,  located  in  an  interval  dx  at  x  and  within  an  angular  interval  d9  at  9. 

If  we  assume  that  the  test  particles  suffer  energy  loss  more  slowly  than 
angular  deflection,  we  can  derive  a  differential  equation  for  f(x, 9) .  If  v  is 
the  same  for  all  particles  then  we  can  write  a  continuity  equation.  The 
divergence  of  f(x,  9)  v  cos9  in  an  interval  dx  di^is  equal  to  the  total  number 
of  particles  scattered  from  all  other  dx’  and  dfl,'  minus  those  scattered  from 
dx  d/1  into  all  dx'  and  d/1'.  This  gives 


(27) 


-  1  f(x, 9)  2v  sin  n  (9')  d/1  '  2Tr  sin9  d9 

Jn'  ® 

where  A'  is  the  solid  angle  of  integration  which  includes  all  scattering 
directions,  which  is  expressed  in  terms  of  .  Since  the  first  integration 
is  to  be  carried  out  over  all  possible  pairs  of  dA  and  dA  it  is  possible 
to  define  a  new  set  of  coordinates  in  which  cylindrical  symmetry  exists.  Thus 
a  factor  of  2ir  can  be  integrated  out  of  the  first  integral.  The  second  depends 
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only  onJl'  and  thus  can  be  expressed  in  terms  of  O',  the  total  cross  section, 
obtain: 

d 


dx 


2n  c^sin  0/2 

'Vos  9 - 


^TTn, 


COS0 


r 


V')o^V')  sin|^  sini^'d/^ 


(28) 


The  integral  term  represents  the  back  scatter;  for  if  we  set  the  integral  term 
equal  to  zero,  then  we  obtain  the  "first  scatter"  solution  fi(x,0) 


fi(x,0)  -  fio  exp 


2s in  9/2 


(29) 


cos  9  .1 

Here  fxo  is  the  number  of  ions  deflected  into  an  element  of  length  dx  at  an 
angle  9;  i .e. 

fio  “  lo  /v  COS0  (30) 

which  gives 

_  Q/o  r  2n„0'sin9/2  I 

-  X  (31) 


.  .  B  COS  9/2 

fi(x,9)  -  Tm^po  ^sirriT;^^ 


COS  9 


which  agrees  with  our  previous  result  Eq.  (I6)  but  includes  the  substitutions 
of  Eqs.  (17)  and  (30) . 

Equation  (29)  can  be  used  to  show  the  buildup  in  density  near  the  body 
by  those  test  particles  which  are  deflected  through  an  angle  near  90’’.  Since 


f(x,9)  -  dn^/d0,  and  in  the  "first  scatter"  approximation 
dn  ■  f  i(  x,  0)  dx 


(52) 


If  we  substitute  9  ■  90°  *  6  then  we  can  write  the  integral  from  0  to  ^9  as 

lAe  ./Fcrne  x/6 

IS  27m  poi*  \  e 


6 


(55) 


where  A9  is  small.  The  choice  of  x  depends  on  the  value  of  A9;  for  x  we 
must  choose  the  minimum  distance  thet  includes  all  the  deflected  particles 
in  the  range  A9  >'6^0  that  have  not  yet  been  deflected.  Therefore  we 


We 
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must  take  for  x 

X  “  AQl/^<rn^ 

If  we  make  the  change  of  variable  z  =  /l©/5  we  obtain 


A6 


r 


CP 


STTn^Po* 


-z  dz 


(5^) 


(35) 


Using  the  Tables  of  Functions  by  Jahnke  and  Emde  we  find 

I  A© 


2Tm  oo®  (0.2194) 


(36) 


To  obtain  a  more  general  solution  for  the  diffusion  equation  we  write 
f(x,  ©)  as 

f(x,©)-  X(x)  @  (©)  .  (57) 

where  this  is  justified  by  the  fact  that  Eq.  (28)  then  becomes  separable 

IT 

®((^)<r(^)sin|^  Sin  d  (58) 

o 

The  integration  of  the  left  side  of  this  expression  yields 

X(x)-Xoe‘''‘  (59) 

This  shows  that  f(x,©)  retains  the  same  dependence  on  x  that  was  found  in  Eq. 
(29).  It  should  be  noted  that  it  is  possible  for  the  coefficient 

of  X  in  the  exponential  factor  to  be  positive.  Thus,  test  particles,  initially 
deflected  forward  by  the  charged  body  may  later  be  scatter  back  so  as  to 
contribute  to  the  density  of  a  certain  region  a  second  time.  The  effect  of 
this  on  the  density,  n,  will  be  considerable  if  po>^  Po-  ^o»  Ps  successive 
scatters  will  remove  the  test  particles  from  the  region  of  compression.  Since 
the  ions  also  slow  down,  thereby,  reducing  the  value  of  po,  the  back  scatter 
is  important  for  oo»  Po  also. 


'  dXLsI. 

X(  x)  dx 


2n^cf'Si 

e 


n  ©/2 


4701* 


cos  © 


@  ( 0) cos© 
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The  function  ®(©)  can  be  expressed  in  terms  of  qo>  the  90®  impact 
parameter  for  encounters  between  field  particles  and  test  particles,  by  the 
express  ion: 


r  r- 


g)  (  0)  ■  ®xp 


J 


2Jrn^qQCos  0/2  +  k  sin®  0/2  sin0  -4n ^sin®  0/2  cos  0/2 


sin*  0/2  (k  COS0  +  2n  ff  sin  0/2) 


-] 


d0 


(40) 


The  Calculation  of  the  Density  as  a  Function  of  Position  in  Cylindrical 
Coordinates 

Let  us  divide  fi(x,0)  by  the  area  into  which  the  particles  at  x  with  a 
deflection  0  will  move.  Since  an  ion  with  impact  parameter  p  will  describe 
a  circle  of  radius  p  that  has  its  center  at  r  where  jl  is  given  by  (see 

Fig.  3)  . 


we  have  for  dn(r,x,0)^0 

dn  .  ft(x.0), 

d0  Tr[(i'  +  p)^  - 


47rp  ypi®^  p®' 


(41) 


(42) 


where  we  have  assumed  that  the  ions  deflected  into  an  angle  0  distribute 
uniformly  over  the  ring  p  <  r  ^  +  p  . 

To  obtain  the  density  at  r,  x  we  integrate  over  all  0  contributing  to 
the  density  at  r: 

n  ®b 

(45) 


n(r,x)  - 


nQ(x,  r,  0)d0 


'0. 


Since  those  ions  scattered  at  p(0  )  will  distribute  out  as  far  as/^(0  )-p(0  ) 

9  d  d 

and  since  this  must  equal  r  if  0  is  to  be  the  minimum  angle  to  contribute 

fl 

to  the  density  at  r  (along  the  positive  x  axis:  i.e.  for  forward  deflections) 
and  since  for  the  largest  angle  to  contribute  0^^  must  be  such  that  i’(9|^)^p(0|^) 

is  equal  to  r.  Therefore,  the  conditions  for  0  and  0.  are: 

A  D 
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0,;  r  -  i(0  J  +  p(9j 

8  8  a 

0k;  r  =  i(ej  -  p{ej 


Substituting  for  fi(x,0)  we  get 


n(x,  r) 


e.  p  fs  sin  0/2  nefi'  1 

b  ngPo  L  cos  0  V _ 

4po  s in0cos0s in®^  ^  Po^cot®  ^  +  po®  sin®  0] 


d0  (46) 


0  ;  r  -  Po®  cot®  0  /2  +  Po®  sin®0 


+  Po  s i n0_ 


®b’  ^  P°^  ^''’^^bj  -  Po  sin  0|^  (48) 

Now  let  us  consider  the  distribution  that  results  if  we  include  the  fact 

'that  the  ions  do  not  distribute  uniformly  over  the  ring  at  x,  r.  Let  us  first 

assume  that  all  the  particles  fi(x,0)  are  deflected  into  one  orbital  circle, 

the  projection  of  the  helix  onto  a  plane  perpendicular  to  x.  The  line  density 

of  particles  is  fi(x,0)/27rp.  Figure  4  shows  that  we  can  write  the  number  of 

ions  in  a  ring  of  thickness  dr  at  r  as 

2  fi(x,0j  _ dt__ 

2irp  |cos  0(1 

where  the  factor  of  2  comes  from  the  two  limbs  of  the  circle  crossed  by  the 
ring.  Now  these  particles  will  actually  be  distributed  randomly  over  this 
element  of  the  ring  (i.e.  the  ring  of  width  dr  at  r) .  Therefore,  the  density 
is  given  by  fi(x,  0)/277®pr  |  cos  0(  |  between  the  points>^-p  and /+p  for  angle  0. 

To  obtain  /cos  a/  we  write,  the  equation  for  a  circle  of  radius  p  that  is 
displaced  from  the  origin  by  a  distance  in  polar  coordinates: 

r®-  2ir  COS0  +  i®  -  p®  (49) 


tan  ® 


therefore 


tan  cx 


/cos0  -  r 
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or 


a=arctan  [cosi-j-^]  ^  ^2^ 

Substituting  for  0  from  Eq.  (^9)  and  from  Eq.  (4l)  we  obtain  an  expression 


for  Qt  which  can  be  used  in  the  density  equation. 

dn(x.  r.9)  fi(x.e) 

dd  2n*prTcosW 


(53) 


If  this  is  integrated  over  9  from  9  to  9.,  where  these  values  are  given  by 

3  D 

Eqs.  (47)  and  (48)  we  obtain: 

re, 


n(x,  r)  - 


2  0/  ^ 

n^Po  cos  /gcxp  j 

sin  ®/2  ne(^ 
cos  9 

1  d9 

g 

27r/oS  i  n9cos9s  i  n®2  rj 

Icos  arctarj 

f  2(po*cot 

‘'/a  -  r*  i| 

.[2r®(po*cot*  ^2po*s  in*9) -po'*cot* 

(5^) 

where  9  and  9.  are  obtained  from  the  solutions  of  expressions  for  r: 
a  b 

e^j  r  ■  I"  po*  cot*  +  po*  sin*  ej  +  p©  sin  9^  (55) 

®b'  ®bj  *  -  po  sin  9|^  (56) 

The  density  given  by  Eq.  (54)  assumes  that  the  plasma  is  cold,  that  the 

velocity  of  the  test  particles  is  constant  and  that  the  test  particles  are 
removed  by  deflections.  The  expression  gives  a  minimum  value  for  the  density 
for  the  case  of  motion  parallel  to  the  magnetic  field. 

6.  Application  of  the  Results  to  the  Calculation  of  the  Compression  of  Ions 
by  a  Satel I i te 

Using  Eqs.  (25)  and  (26)  we  can  obtain  the  value  of  the  average  density 
of  the  deflected  ions.  Using  9  •  90°  end  9  ■  120®  in  Eq.  (26)  and  setting 
n  equal  to  10®  lons/cm®,  then  we  get,  including  the  factor  2  for  back  scatter, 
a  value  of  n  ■  4.0  x  10®  lons/cm®  for  H*  ions.  For  0*  ions  the  value  is, 
however,  small.  The  value  of  n  must  be  reduced  somewhat  since  p  is  not 
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sufficiently  small.  If  this  is  taken  into  account  and  if  we  include  the  back 
scatter,  then  we  obtain  foi  n  for  ions  in  ions/cm^ 


a  =  1  m  8  X  10®  3.2  X  10®  4.0  x  10® 


Since  the  compression  effect  is  negligible  for  0  ions  but  large  for 
ions,  satellite  compression  effects  will  be  masked  below  altitudes  of 
about  600  km.  Since  below  this  altitude  the  density  of  0*  ions  is  10  times 

4*  4* 

as  great  as  the  density  of  H  .  The  density  of  the  0  drops  off  much  more 
rapidly  than  the  so  that  their  densities  are  equal  at  an  altitude  of  ,1,500 
km  and  dominates  above  that  altitude.  Therefore,  large  satellite  wakes 
would  be  masked  below  about  6OO  km. 

The  scale  length  of  the  density  distribution  will  be  x©  which  is  the 
mean  free  path  \  times  the  facto .  Using  an  average  value  of 
8  ■  120®,  this  factor  becomes  O.3I8.  From  Table  I  we  see  that  for  H^, 
xq  ■  1 .6  km  and  for  0^,  x©  ■  290  km. 

7.  Deflection  by  a  Charged  Sphere  Passing  Through  an  Ionized  Gas  Cloud 

So  far  we  have  assumed  the  charged  body  to  be  moving  through  an  ionized 
gas  of  uniform  density  and  infinite  extension.  With  these  assumptions  we 
obtained  a  solution  for  fj. 

Let  us  now  assume  that  the  charged  body  moves  relative  to  a  distribution 
of  ionized  gas  of  density  t  +  x/v)  where  x  =  0  is  the  position  of  the  body 
at  t  ■  0,  and  we  are  in  the  system  of  coordinates  moving  with  the  charged  body. 
Therefore  we  obtain  from  Eq.  (28),  neglecting  the  integral  term, 

r  2  sin  I  ~l 

fi(x,©,t)  ■  fxo  exp - ~  Q  \  njt  +  x/v)dxl  (57) 

The  function  fi©  is  exactly  the  same  as  before  except  that  we  must  substitute 
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for  the  function  n^(  t  -  x/v  cos  9).  The  resulting  equation  is,  therefore 


fl(x,9,t)  -os0°si'n^\/2  “ 


V  COS0 


exp 


f  ^  1 

2s  in  0/2 

\  ng(  t  +  X/ v)  dx  / 

COS  6 

L 

(58) 


0 


8.  Derivation  of  the  Density  of  Ions  for  the  Motion  of  a  Satellite  Through 

a  Cold  Plasma  at  an  Angle  to  the  Magnetic  Field. 

We  consider  the  case  where  the  velocity  vector  of  the  satellite  makes 
an  angle/3  to  the  magnetic  field.  We  have  already  obtained  an  expression  for 
f(x,0)  for  the  first  scatter  approximation.  We  will  drop  the  exponential 
factor  in  the  present  derivation  thus  neglecting  the  scatter  of  test  particles 
by  the  field  particles  and  consider  only  fio(x,S‘)  where 

fio(x,e')  l-(^s0fTinaVv2  (59) 

This  also  neglects  back  scatter,  however,  since  \  is  large,  back  scatter  will 
be  unimportant  when  the  satellite  moves  at  an  angle  to  the  magnetic  field. 

For  the  present  problem,  we  must  take  x  parallel  to  B.  Thus,  the  cos9 

term  will  remain  unchanged  in  the  presence  of  a  magnetic  field  since  it 

arises  from  the  value  of  the  component  of  motion  of  the  test  particle 
parallel  to  the  magnetic  field.  However,  the  angle  0‘  is  the  deflection 
angle.  In  order  to  transform  fio(x,9' ) -  into  a  form  appropriate  to  the  new 
orientation  of  the  magnetic  field,  we  will  introduce  the  azmuthal  angle  0'. 
Define 


G,n(x  9'  0')  -  .  nS-SsL.  ,  cos  9'/.2 

Ql0^x,9  2ICOS0I  sin=*0'/2 


(60) 


We  now  need  to  transform  Gio(x,9',0')  into  Gxo(x,9",0")  where  d"  and  0"  are 
the  coordinates  in  the  lab  system  (see  Fig.  ^) .  Since  the  velocity  of  the 
coordinate  system  is  the  same  as  that  of  the  deflected  test  particles  we  have 
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that 


Thus 


011  „  90°  +  9' 


sin  9“ 


(61) 

(62) 


Gxo(x,e")0")  “  2  lcos9|  cos=*  0" 

Here  again  the  |  cos  @|  term  remains  unaltered. 

We  will  now  transform  from  the  coordinates  x,9",0"  measured  relative  to 
the  axis  of  motion  to  z,0,0  measured  relative  to  the  direction  of  the  magnetic 
field  vector  (see  Fig.  6).  The  factor  sin9"  cos  ^0"  must  now  be  transferred 
into  the  coordinate  system  0,0  as  measured  about  Z.  Using  the  law  of  cosins 
for  a  spherical  triangle  we  obtain 


s in  0"  ^  II  -  ( COS0  cosP  sin0  sinP  cos0)^J 

cos^0"  "  icos0  cosP  +  sin©  sin3  cos0j  ® 


(65) 


We  obtain  for  Gio: 


GUz  9  a)  •  1  t  ~  5in0  Sinp  ,cps0}«7  ‘ 

'  '  2  |cos6f  [cos©  COS0  +  sin0  sinp  cos^J'*  ' 


(65) 


(66) 


The  velocity  of  the  test  particle  in  the  lab  system,  v',  is 
v'  ■  2v(cos0  cosp  +  sin©  sin0  co50) 

Consider  now  the  distance  st  over  which  the  charge  distribution  will  be  dis¬ 
tributed  in  a  time  interval  t.  From  Fig.  7  we  have 

s^  =  V^  [^1  +  4cos^0(  cos©  cosp  +  sin©  sing  cos0)® 

-  41cos0lcosj3  (cos0  cosg  +  sin0  sing  cos0)J 
Now  the  test  particles  that  move  parallel  to  6  with  a  forward  velocity 
v' |cos0|  will  be  spread  out  along  the  line  st.  Therefore,  the  function 
Gio(z,0,0)  must  be  multiplied  by  the  factor  v'  |  cos@//s.  Thus 

g!io(  l5)  •  2(  cos0cosg+s  1  n0s  i  ngco50/  [1  +4(  cosOcosg+s  i  n0s  i  ngcos0)®cos®9 
-  4(cos9cosg+sin0singcos0) |cos0|  cosgj  (67) 
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Now  to  obtain  the  density  we  divide  fio(z,6)  by  v(cos0,  -  cosSq) t  and  integrate 
over  A®  (i.e.  we  divide  by  vt  sin6)  and  also  divide  by  the  width  of  the  distri¬ 
bution  which  is 

^p(  9)  =  ^Po  s  i  n9  (  68) 

where  p»  p  or  by 

2p(e')  =■  2po  cot©"/2  (69) 

for  p>’5‘  p-  Thus  we  get 


f  ( z,  9)  = 

10 


2Tr 


9io(  z,  9,  l3)  d(5 


(70) 


so 


,  Q  _ _ fi.r(Z|  6} _ 

'  ’  ’  ^p(9  )  vt  sin© 


for  p»  p  and 


i/a 


7  a  rl  s  _ ^1iQ( 

nu^e,  t;  2p(e')  vt  sin® 
for  p>^  p.  If  we  transform  p(0')  we  get 

,  ,  [l  -  (cos9  cosP  +  sin®  sinP  cosl8)^-7 

(cos9  cosP  +  sin9  sinP  cos0) 

The  results  of  the  substitutions  of  Eq.  (67),  (68),  (70),  and  (75)  into  Eq. 

(71)  and  (72)  give  for  p»  p 


(7i) 


(72) 


(75) 


nfx  ©  tl  >  _ 

4po  Sin*  9  vt 


(74) 


and  for  p  p 


n(x,  9,  t) 


fin(x.9)(cos9  cosP  •«•  sin®  sin3  cost)) 
2po  vt  [^l-(cos0  cosP  +  sin®  sinP  cosl5)*J 


(75) 


where  ^ 

i'io(z>®)  ■  I  Tco^l  *  coslS)^  ^^|"(cos0  cosP 

■^0 

+  sin®  sinP  cos£))*  fl+4(co50  cosP  +  sln9  sinP  cos0)*  cos*  9 
-4(cos®  cosP  +  sin®  sinP  cos0)  |  cos0|  cosfij  J  d0  (76) 


I 
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The  distribution  that  is  produced  gives  rise  to  a  ^shock  front*!  If  the 
angles  between  these  fronts  and  the  velocity  vector  of  the  charged  body  are 
1  and  72>  (see  Fig.  2)^  then  these  angles  are  given  by  the  expression 


Let  us  now  determine  the  angle  for  which  n(Z,9,t)  will  be  the  largest, 
and  then  determine  how  far  the  high  density  region  extends  away  from  the  body 
before  the  density  of  test  particles  plus  the  density  of  field  particles 
becomes  2n^.  For  p>>  p,  Eq.  (75)^  there  is  an  infinity  produced  for 

1  -  ( COS0  cosP  +  sin0  sing  cos0)®  ■  0  (80) 

However,  this  is  spurious  since  it  arises  from  the  breakdown  of  the  approximation 
p>>  p.  Likewise,  the  infinity  produced  by  9  =  90°  is  spurious,  arising  from 
the  neglect  of  the  exponential  decay  function. 

The  integration  over  0  presents  considerable  difficulty  when  compared 
with  the  physical  value  of  such  an  integral  It  is  sufficient  to  note  that 
the  integration  will  produce  a  slight  smoothing  out  of  the  shock  front.  We 
should  note,  also,  that  the  factor  f io{Z,9} /2p  may  be  replaced  by  n^/po  where 
is  the  density  of  test  particles.  Since,  for  a  positive  sphere,  values 
for  7  of  about  10  can  be  obtained,  we  have 


°a:£o  .  n 

vt  e  vt 


(81) 
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Thus  the  density  of  test  particles  will  equal  the  density  of  fieid  particles 
when  vt'^^lOpo*  If  Po  's  3  meters,  then  this  gives  a  wake  that  extends  JO 
meters  before  the  density  has  dropped  to  twice  the  ambient.  If  v II  B  then 
for  ions  and  an  impact  parameter  of  3  meters,  the  cross  section  of  the 
wake  would  be  about  5  x  10'*  m^.  Because  of  the  geometry  of  the  wake,  a  long 
rod  shaped  region,  the  effective  radio  cross  section  would  be  larger. 

The  reflection  of  a  radio  signal  can  be  either  coherent  or  non-coherent . 
If  non-coherent,  then  we  would  need  an  impossibly  large  number  of  electrons, 
more  than  10^°.^  If  the  reflection  is  coherent  then  it  can  be  either 
metallic  or  non-metal  I i c.  Since  the  index  of  refraction  for  a  wave,  n*, 
of  frequency  f,  is 


n 


*  S 


(82) 


(in  cgs  units)  we  obtain  metallic  reflections  >4ien 

I0“e  (95) 

Thus,  when  f  is  2  x  10^  sec  ^  we  require  n  ■  4  x  10®.  Thus  we  should  be  able 
to  obtain  large  radar  cross  sections.  For  the  above  example  the  cross  section 
is  about  180  m®,  including  both  halves  of  the  shock  front. 

9.  Deflection  of  Ions  and  Electrons  by  a  Screened  Electric  Field. 

Thus  far  we  have  considered  the  charged  body  to  be  unscreened.  This  is 
accurate  either  for  charged  bodies  in  a  low  density  plasma  or  for  very  small 
bodies.  For  satellites  such  approximations  are  not  satisfactory.  The 
screening  for  probes  an  for  satellites  has  been  extensively  considered  by 
Walker*.  The  potential  falls  off  very  rapidly  down  to  a  value  of  about  kT; 
however,  the  drop  is  several  times  slower  than  that  given  by  the  assumption 
of  local  thermodynamic  equilibrium.  Below  kT  the  field  falls  off  more  slowly 
than  r  .  Even  for  larger  values  of  the  impact  parameter  for  which  the 
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potential  will  be  only  about  .1  kT  significant  compression  of  the  test  particles 
can  occur. 

There  will  also  exist  in  this  screened  electric  field  a  pericritical 
surface  for  which  particle  orbits  become  critical  and  the  particle  spirals 
into  the  charged  body.  We  must  also  consider  the  fact  that  most  all  of  the 
particles  that  strike  the  surface  of  the  satellite  will  be  neutralized.  Since, 
for  a  negative  sphere  all  of  those  ions  that  reach  the  pericritical  surface 
will  be  accreted  to  the  surface  offthe  charged  body,  the  ions  with  small 
impact  parameters  cannot  be  important  in  the  production  of  large  wakes.  For 
a  positively  charged  body,  the  ions  will  be  important  in  the  production  of 
large  wakes  but  the  electrons  that  reach  the  pericritical  surface  will  be 
accreted  and,  for  the  most  part  neutralized  at  the  surface  of  the  body.  For 
small  charged  bodies,  such  as  antennae  moving  at  satellite  velocities  the 
pericritical  surface  may  become  distorted  due  to  the  relative  motion  of  the 
ions.  Thus  the  focal  point  of  the  spiral  orbits  will  be  shifted,  and  many 
of  the  ions  will  not  be  accreted.  It  must  be  noted,  however,  that  in  the 
presence  of  a  magnetic  field,  the  screening  process  may  be  sufficiently 
altered  to  allow  a  greater  number  of  particles  to  be  compressed  without 
their  being  accreted  to  the  body. 

Thus,  in  most  cases,  it  will  be  the  particles  with  the  same  sign  as  the 
satellite  which  will  be  compressed.  The  result  will  be  that  a  cloud  of  charge 
will  be  produced  in  the  neighborhood  of  the  body,  and  will  move  away  along 
the  magnetic  field  lines.  If  the  body  moves  parallel  to  the  magnetic  field 
an  extensive  cloud  of  charge  will  be  produced  that  moves  ahead  of  and  behind 
the  satel 1 i te. 

Using  these  facts,  we  can  obtain  descriptions  for  the  form  of  the  wake 
under  various  circumstances.  If  a  negatively  charged  sphere  moves  with 
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satelllte  velocities  parallel  to  the  magnetic  field,  the  ions  will  be  accreted 
to  the  surface  of  the  body  and  will  not  be  compressed.  The  electrons,  because 
of  their  velocity,  will  be  reflected  on  both  the  front  and  rear  sides  of  the 
body.  However,  since  the  electrons  that  are  reflected  merely  take  the 
place  of  those  that  normally  would  come  from  the  opposite  side  of  the  body, 
there  will  arise  only  a  slight  compression  due  to  the  displacement  of  the 
particles.  If  the  body  were  moving  with  a  velocity  comparable  to  the  electron 
thermal  velocity,  the  compression  of  the  electrons  would  be  significant.  If 
the  negatively  charged  sphere  is  moving  at  an  angle  to  the  magnetic  field, 
essentially  the  same  thing  will  occur.  The  Ions  are  accreted  to  the  surface 
of  the  body  and  neutralized  while  the  electrons  are  reflected  and  produce  only 
a  si  Ight  compression. 

For  «  positively  charged  body  moving  parallel  or  at  an  angle  to  the 
magnetic  field,  the  electrons  will  be  accreted  to  the  surface  of  the  body 
while  the  Ions  are  deflected  and  produce  a  wake  with  a  density  several  timas 
the  ambient  density. 

Still  other  factors  may  be  Important!  If  tha  body  Is  equipped  with  Ion 
collectors  or  other  apparatus  that  Is  driven  positive  with  respect  to  the  rest 
of  the  body,  then  this  portion  of  the  vehicle  will  produce  plasma  compressions 
In  the  manner  described. 
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0  ions 

ions 

Thermal  Velocities 

1 .55  km/sec  . 

Satel 1 itc  Velocities 

7  km/ sec 

Thermal  Velocities 
6.1  km/sec 

Satel iite  Velocities 

7  km/ sec 

Po 

25.6 

1  .24 

25.8 

19-5 

Po 

8.52 

39-0 

2.12 

2.44 

m 

9x10® 

5x10® 

h 

8.1+5  X  lO"® 

Table  I.  Important  parameters 

The  values  of  po^  P<,i  h  singly  ionized  oxygen  and  hydrogen 
(  Z  ■  I)  :  Satel  I  1  te  potent  ial  V  >  10  vol  ts,  satel  I  1  te  radius  a  I  m, 
ambient  electron  density  ■  lO*-*’  m"®,  temperature  T  ■  1500°  K,  magnetic 
induction  B  ■  5  x  10"°  w/m®,  A  is  the  atomic  weight,  v  is  the  velocity  of 
the  ion  and  Vj_  is  the  component  of  the  velocity  perpendicular  to  the  magnetic 
field.  The  results  are  expressed  in  meters. 
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Fig.  1. 

Fig.  2. 

Fig.  3. 

Fig.  4. 

Fig.  5. 

Fig.  6. 

Fig.  7. 


Diagram  showing  parameters  entering  into  the  trajectory  of  a  cold 
proton  moving  in  the  electric  field  of  a  charged  body  as  seen  in  the 
center  of  mass  system. 

The  form  of  the  shock  fronts  formed  by  the  passage  of  a  charged  body 
through  a  cold  plasma  at  an  angle  to  the  magnetic  field. 

Diagram  showing  the  relation  of  to  p,  the  radius  of  gyration  and  p, 
the  impact  parameter. 

Diagram  showing  parameters  entering  into  the  calculation  of  the  density 
of  test  particles  in  a  ring  of  thickness  dr  and  radius  r. 

Illustration  of  the  relation  between  the  lab  coordinates  and  the  center 
of  mass  coordinates. 

Transformation  of  the  coordinates  x,  ©*',  0"  measured  relative  to  the 
direction  of  the  magnetic  field  vector. 

Illustration  of  the  distance  over  which  the  charge  cloud  is  distributed 
by  the  motion  from  the  point  of  initial  deflection  to  a  position 
v' /cos9/t. 
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REFLECTION  OF  RADIO  WAVES  INCIDENT  ON  BOUNDARY  SURFACES  WITH  APPLICATION  TO 


THE  THEORY  OF  THE  KRAUS  EFFECT 
E.  H.  Walker 

Abstract 

Satellite  wakes  present  sharp  boundaries  to  radio  waves  that  are  incident 
upon  them.  Also  the  wake  produced  by  a  satellite  is  not  free  to  dispurse  in  all 
directions  but  is  constrained  by  the  earth's  magnetic  field  so  as  to  form  a  large 
and  sharply  defined  boundary  surface  which  persists  long  after  the  passage  of  the 
satellite.  In  order  to  show  the  relation  between  such  satellite  wakes  and  the 
Kraus  effect  a  treatment  of  the  reflection  of  radio  waves  off  sharp  boundary 
surfaces  is  presented.  Although  these  wakes  have  been  considered  by  many  authors 
to  be  insufficient  to  account  for  the  Kraus  effect,  we  show  that  the  reflection 
coefficient  for  these  boundaries  become  significant  under  favorable  conditions 
of  obstrvation.  The  very  large  size  of  the  wakes  and  their  geometry  result  In 
an  inverse  square  law  rather  than  an  inverse  fourth  power  law  relating  the 
power  to  signal  strengths  for  radar  observations. 


REFLECTION  OF  RADIO  WAVES  INCIDENT  ON  BOUNDARY  SURFACE  WITH  APPLICATION 
TO  THE  KRAUS  EFFECT 

1.  Introduction 

In  a  previous  paper^  we  presented  a  mechanism  involving  the  deflection 
of  the  ions  of  the  plasma  by  a  charged  satellite  so  as  to  form  a  wake.  The 
earth's  magnetic  field  wiil  produce  a  compression  of  this  plasma  in  the 
immediate  neighborhood  and^  If  the  satellite  moves  parallel  to  the  magnetic 
field  an  extensive  region  of  compressed  plasma  will  result  in  favorable  cases. 
Because  of  the  constraint  of  the  earth's  magnetic  field  and  because  of  the 
large  mean  free  path  of  the  ions  in  the  upper  Ionosphere,  these  ions  maintain 
a  coherent  geometry  long  after  their  density  has  dropped  to  very  low  values, 
thus  presenting  large  sharply  defined  boundary  surfaces. 

It  is  most  important  that  these  surfaces  are  sharp  and  flat  (or  at 
least  with  only  a  slight  curvature).  Due  to  random  fluctuations  in  density, 
regions  of  high  or  low  density  constantly  appear  throughout  the  Ionosphere. 
However,  those  regions  are  of  small  size  and  occur  with  a  random  distribution 
and  with  various  random  shapes.  Thus,  even  though  the  wake  of  the  satellite 
is  superimposed  on  a  medium  having  fluctuations  in  density,  the  fact  that  the 
reflection  from  the  wake  is  coherent  whereas  the  reflection  from  the  random 
fluctuations  is  incoherent  causes  the  satellite  related  radio  reflections  to 
be  much  stronger,  when  they  occur,  than  in  background  reflections  from  the 
random  density  fluctuations. 

This  distinction  between  coherent  and  incoherent  reflections  is  obvious 
from  a  consideration  of  the  expression  for  radiation  from  an  accelerated 
charge.  If  a  point  charge  q  undergoes  an  acceleration  u,  then  the  rate  at 

which  this  body  radiates  energy  in  the  form  of  electromagnetic  radiation, 

a 

dW/dt  is 


■2 


dW  cr  u  / ,  < 

dt  “  Co  ^  ^ 

The  dependence  upon  the  total  amount  of  charge  Is  very  strong,  varying  as  the 

square  of  the  total  charge  that  Is  radiating  coherently.  For  a  collection  of 

charges  that  are  distributed  over  a  region  of  space  and  which  radiate  out  of 

phase  with  one  another,  the  total  radiation  will  be  greatly  reduced,  no  longer 

depending  upon  the  square  of  the  total  charge.  The  amount  of  radiation  from 

such  an  assemblage  of  charges  depends  upon  their  geometry  and  phase  relation 

but  these  might  be  so  chosen  as  to  produce  only  quadrapole  or  a  higher 

multipole  radiation. 

When  the  terms  “coherent'*  and  “incoherent"  are  applied  to  the  reflection 
and  scattering  of  electromagnetic  waves  from  a  region,  they  indicate  whether 
the  wave  was  reflected  by  charges  that  are  reradiating  in  phase,  or  by  charges 
that  are  out  of  phase. 

The  following  sections  give  a  treatment  of  the  reflection  of  radio  waves 
off  of  sharp  boundary  surfaces  that  separate  regions  of  plasma  that  are  at 
different  densities.  We  will  first  consider  the  problem  quite  generally  and 
then  apply  this  to  the  Kraus  effect. 

2.  Derivation  of  the  General  Reflection  Equations. 

Let  us  first  consider  the  motion  of  an  electron  in  a  plasma  when  it  is 
subjected  to  an  electric  field  E. 

The  drift  velocity  Vq  of  an  electron  with  mass  m  and  charge  e  that 
moves  freely  except  for  collisions  with  ions  and  other  electrons  can  be 
related  to  the  electric  field  by  the  equation 

''o  (2) 


where is  the  average  time  between  colissions.  Since  the  current  density 


j  ■  ne  V. 


(3) 


where  n  is  the  density  of  electrons,  we  have  the  relation  between  E  and  j: 

ne^E  _  di  .  1  . 
m  dt  T 

The  Maxwell's  equations  can  be  used  to  obtain  the  equation 


^  E  •M 


St 


If  we  now  substitute  into  Eqs.  (k)  and  (^) 

•icut 


we  obtain 


and 


J  “  Joe 


ne^  E  I  "  luiT  . 


o  b^r 

V*E  »  ~  \uiMj  *■  n€  gp- 

Thus,  upon  substituting  Eq.  (7)  into  Eq.  (8)  we  have 


.  2-  ia)-r  pne*  ,  ^ 

V  E  ■  '  i"  ^ — •  E  +  [>(. 

^  I  “  iu)t  m 

If  we  substitute  for  E  the  expression 

E(l)  =  Eo(6) 
into  Eq.  (9)  we  have 


8=E 

5? 


iU)T  UfiC: 


t  ■  peu  2  E 


1  -  iWT  !ll 

In  metals  tw  I  for  frequencies  below  the  optical  region.  This 
mation  along  with  the  relation  between  t  and  the  conductivity  c: 


0  = 


_  QSr.T 


can  be  used  in  Eq.  (II)  to  obtain 
gpr  s  -  p«w*E  -  iwpaE 

This  last  equation  governs  the  so-called  metallic  reflection  phenomena. 


(^) 

(5) 

(6) 

(7) 

(0) 

(9) 

(10) 

(11) 

approxi- 

(12) 

(13) 
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If  we  assume  -no  »  I  then  Eq.  (II)  becomes 

E  ( lU) 

This  expression  applies  In  those  cases  where  the  electrons  can  be  considered  as 
non-interacting  with  other  particles.  We  will  refer  to  this  type  of  reflection 
as  non-metal  1 ic. 

In  the  ionosphere,  both  metallic  and  non-mettalic  reflections  can  occur. 
The  angular  frequency  <>>  that  divides  the  two  types  of  reflections  depends  upon 
T.  According  to  Spitzer^,  the  resistivity  p  is 

P  “  65.3  (ohm  -  m)  (I5) 

where  T  is  the  absolute  temperature  and  the  value  JLnJV  can  be  obtained  from 
Table  I .  From  Eq.  ( 12) 


Values  ofi>7yLfor  various  values  of  T  and  n  according  to  Spltzer^.' 


We  obtain  for  t 


r»/a 


•  fer  •  io»  ..c 


(16) 


where  n  is  the  number  of  electrons  per  cubic  meter. 
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In  a  region  of  the  ionosphere  where  the  density  is  m  ^  and  the  temperature 

is  1000  we  obtain  for  the  conditions  on  u  for  metallic  reflection 

w  «  7.»t5  •  10=  sec  ^ 

i.e.j  for  metallic  reflection  to  dominate  the  frequency  of  the  radio  wave  must 
be  much  less  than  Il8  cycley^ec  the  reflections  are  non-metallic.  The  penetra¬ 
tion  depths  for  both  metallic  and  non-metallic  reflection  from  the  ionosphere  are 
small  compared  to  the  thickness  of  the  ionosphere  as  will  be  shown  in  section  4. 

Coherent  and  Incoherent  metallic  reflections. 

Equation  (13)  can  now  be  used  along  with  the  boundary  conditions  to  obtain 
a  complete  description  of  metallic  reflections.  The  boundary  conditions  that 
relate  E  of  the  Incident  ray,  E“  of  the  reflected  ray  and  E*  of  the  transmitted 
ray  are 


£  X  (£  +  £")-  n  X  i* 


(17) 


and 

J1  X  (Js  X  i  +  X  £")/Mi  -  n  X  (V  X  i')/HB  .  (18) 

where  £  is  the  unit  vector  normal  to  the  surface  and  £  is  the  propagation  vector. 
Similar  expressions  relate  the  magnetic  field  vectors  H,  H*'  and  H.'  across  a 
boundary. 

If  we  express  Eq.  (13)  as 


=  -  (iA€w=  +  i(i»|i>r)E  =  -  K=E 


where 


K .  a .  le .  k  /TTE 


Then  the  boundary  conditions  will  give 


and 


£--£o'  ei(“S-ut) 


u.  .  £  ^  £'  C0».,+  IP.) 

“  MW  k 


(19) 

(20) 

(21) 


(22) 
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Furthermore,  the  reflection  coefficient  r  can  be  shown  to  be 


E" 

2 

1  -  02  -  6®  -  2iB 

E 

1  +  aa  +  OP®  +  B* 

(23) 


in  this  general  situation. 

The  condition  that  the  reflection  be  coherent  is  equivaient  to  setting 
eu)  «  o  ^  a  condition  that  Is  met  in  metals  below  the  optical  region.  In  this 
case,  Eq.  (21)  and  (22)  become 


£1  a  (J  I  -  Wt) 


(24) 


and 


H'  =  (1  +  i) 


ijr 


(25) 


These  fields  fail  of  to  i/e  of  their  surface  value  in  a  distance  the  skin 
depth 


6 


I 


(26) 


The  expression  for  the  reflection  coefficient,  Eq.  (23),  under  this  approximation 
becomes 


-  2  (27) 

4,  Non-Hetallic  Reflections 

In  a  low  density  plasma  the  expression  for  metallic  reflection  are  no 
longer  valid  so  that  we  must  return  to  Eq.  (14).  We  will  express  this  in  the 
form 


E  -  -  E  (28) 

where 

K  peu*  -  (29) 

We  thus  have 

i'  -  JEa  (50) 
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(5.) 

If  -  <  3. 17-  10®  (m.K.S.  units)  then  this  wave  is  damped  out  in  a  characteristic 

distance 


•  /T;  a  i 

PS.® 

V  m 

A  calculation  of  the  reflection  coefficient  yields  for  a  wave  going  from  a  region 
of  density  ni  to  a  region  of  deni  sty  ns 


It  Is  interesting  to  note  from  Eqs.  (26)  and  (32)  that  under  the  conditions 

of  the  ionosphere  at  T  =•  1000°K  and  n  =  10^  m  ®  and  at  the  frequency  of 

I18  cycles  sec  that  divides  the  metallic  reflections  from  non-metallic  reflections,  & 
has  a  value  of  7*55  m  while  6'  has  a  value  of  m.  indeed,  the  ratio  6/6' 
always  equal  at  the  frequency  where  ut  «  I  so  long  a  n>  2*10*m  ®. 

3.  Application  to  Radio  Reflections  in  the  Ionosphere 

In  the  region  of  the  ionosphere,  reflections  are  non-metallic  except  at 
very  low  frequencies.  Thus,  the  reflections  of  signals  at  all  but  the  lowest 
frequencies  will  be  governed  by  Eq.  (33)'  Let  us  consider  a  few  typical  cases 
of  reflections  in  the  ionosphere.  If  a  20  Hc^ec  wave  travels  from  empty  space 
into  a  region  of  density  lO^^m"®  (or  from  one  region  at  a  density  of  lO^^m  ®,  m, 

into  one  of  density  2* 10^^  cm~®,  ng),  then  Eq.  (33)  can  be  approximated  as 

('’8  -  ni)*  {'^k) 


thereby  giving  a  value  of 


r  -  2.57  •  10 
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If  the  densities  are  10  times  as  great  then  r  »  207 ' 10 

Let  us  now  consider  the  case  where  S? ^  }  i.e.,  let  us  consider 

the  case  of  reflections  where  n  is  nearly  equal  to  the  critical  frequency.  Thus 
we  can  define  a  quantity  f  such  that 

f  -  1  -  •“►<<  1  (35) 

where  n^  is  the  critical  density  which  is  equal  to  3*15  •  10  *  (m  ®).  We 
will  assume  that  ng  Is  Just  slightly  greater 

r\8  ■  ni  +  <Sn  (56) 

so  that  in/ni  f.  We  obtain: 


As  an  Manple,  take  n  to  be  1^  below  the  critical  density  and  in/ni  ■  lo'®. 

The  reflection  coefficient  for  this  case  is  I/1600. 

This  behavior  is  quite  important  to  the  understanding  of  radio  reflections 
from  satellite  wakes.  As  an  example,  consider  a  satellite  of  radius  I  meter 
passing  a  particular  position  in  the  earth's  magnetic  field  and  producing  locally 
a  compression  of  the  ions  to  twice  their  ambient  density  (one  can  also  consider 
the  results  of  a  satellite  sweeping  out  the  ions  with  the  same  results).  Since 
the  ions  spiral  along  the  magnetic  field  lines,  as  shown  in  our  previous  paper^, 
they  are  allowed  to  spread  out  only  in  one  direction.  Thus,  when  the  ions  have 
traveled  100  km,  they  will  be  confined  within  a  region  of  about  10  km  along 
the  magnetic  field  lines.  Thus,  within  this  region  the  ion  density  would  be 
enhanced  by  an  amount  8n/nx  I0~®.  In  a  region  with  a  density  within  1^  of 
the  critical  density,  r  ■  (l/400)‘.  This  Value  in  Itself  is  small,  but  when 
we  consider  the  other  significant  factors,  such  reflections  become  very  important. 
Since  the  reflecting  surface  is  not  small  compared  to  the  wavelength  of  the 
incident  radiation,  but,  on  the  contrary,  is  quite  large  with  the  shape  of  a  flat 
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sheet  (or  of  a  series  of  rods,  if  the  wake  fluctuates),  a  reflected  radar  signal 
may  fall  off  as  R”*  or  R"®  rather  than  R  *  as  is  usually  the  case.  So  far, 
however,  we  have  not  considered  the  dependence  of  r  upon  the  thickness  of  the 
sheet. 

6.  Non-Hetallic  Reflection  from  a  Layer  of  Electrons 

Consider  a  slab  of  free  electrons  of  thickness/^  with  an  incident  electro¬ 
magnetic  wave  on  the  surface  a.  Reflection  occurs  at  this  surface.  A  second 
reflection  occurs  on  the  second  boundary  surface  b  with  a  portion  being  trans¬ 
mitted  out  of  the  slab.  In  Fig.  I  the  various  parameters  are  explained;  E  is 
the  incident  field,  E"  is  the  total  reflected  field,  E'"  is  the  total  field  trans¬ 
mitted  out  of  the  slab,  and  Ei,  Es  are  the  fields  inside  the  slab.  The  value 
of  Ex  at  the  surface  b  is  related  to  the  value  at  the  surface  a  by  the 


expression 

Ei  (b)  -  Ei  (a)  s  (38) 

From  the  boundary  equations,  Eqs.  (I?)  and  (18),  we  obtain 

E  +  E«  -  Ei  ♦  TEa  (39) 

TEi  +  Ei  =  E"  (40) 

kE  -  kE"  =  k'  El  -  yk'  Ei  (4l) 

Sk'  Ei  -  k'  Ez  =  k  E'"  (42) 

Substitution  of  Eq.  (39)  into  (4l)  to  eliminate  Ex  yields 

kE  -  kE"  -  k'(E  +  E")  -  23'k'  Ei  (43) 

Substitution  of  (39)  into  (40)  to  eliminate  Ex  yields 

y(E  +  E"  -  JEi)  +  Ei  =  E'"  (44) 

Substitution  of  (40)  into  (42)  to  eliminate  Ex  yields 


(45) 
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Substitutlon  of  (4V)  Into  (4^)  to  eliminate  t"'  yields 
;r(E  +  E")  =  [y*  -  1  +  Ez 

Substitution  of  (46)  into  (43)  yields  upon  simpllficati 


E" 

2 

2kk'  -  l<?( 

1  +  y®)-k'®( 

.1  7 

r  * 

E 

X 

2kk'(l  +Jr 

)  +  (k'“  +  k 

on 

2 


(46) 


(47) 


Since  for  our  case  k'  >  K  as  given  in  Eq.  (29)>  by  defining  cX  >  K/k  we  can  write 
Eq>  (47)  more  simply  as 


got  (1  , 

2a  (i  +  /*)  +  0  +ot^)(i  - 


(48) 


If  we  assume  that  ri2  =  nx  +  cfn  and  that  nx  is  sufficiently  close  to  the  value  of 

n^  so  that  the  inequality  of  Eq.  (3^)  holds  then  we  get>  where  K/  is  small: 

2 

(49) 


r  -  (^)  (1  +4?  +  i*) 


fn 


This  expression  yields  essentially  the  same  results  that  were  obtained  in  Eq.  (37) « 
Thus  we  see  that  for  a  thin  layer  of  free  electrons  (where  is  small  but  large 
enough  so  that  the  region  can  be  considered  homogeneous),  the  reflection  coefficient 
is  nearly  Identical  to  the  value  for  an  infinite  layer.  Reflections  by  a  succession 
of  such  layers  would  be  rather  accumulative  if  the  separation  distance  between 
the  layers  allowed  the  successive  wave  fronts  to  be  in  phase. 

Thus  under  proper  conditions  a  satellite  would  produce  a  slab-like  region 
(or  a  series  of  rod-like  regions  if  the  wake  has  an  oscillatory  nature)  of 
slightly  higher  density  than  its  surroundings,  which  are  assumed  to  have  a 
density  that  is  near  the  critical  frequency.  This  region  would  be  highly 
directional  as  a  reflector  of  radio  waves.  The  present  treatment  is  still  too 
general  to  obtain  specific  values  of  the  cross  section.  Nevertheless,  the  values 
should  be  high. 
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7>  Suggestions  For  Observational  Techniques 

These  considerations  at  once  suggest  some  useful  modifications  of  the 
current  experimental  systems.  The  region  from  which  radio  reflections  of  satellite 
wake  phenomena  are  to  be  expected  is  limited  to  the  region  for  which  the  ambient 
density  of  electrons  almost  reaches  the  critical  value.  This  region  should  be 
mapped  out  during  the  experiment  by  means  of  radio  soundings.  The  experiment 
should  be  conducted  only  when  this  region  is  limited  in  extent  and,  furthermore^ 
when  known  objects  will  be  passing  through  the  region.  The  early  morning  hours 
should  be  best  for  this  purpose  since  observation  would  be  made  from  a  region 
of  low  density  toward  a  nearby  region  of  high  density. 

Since  the  orientation  of  the  reflecting  surface  is  determined  by  the 
satellite  trajectory  and  the  magnetic  field;  one  can  choose  the  optimum  geometry 
for  the  radar  transmitter  and  the  several  receivers. 

A  satellite  In  an  orbit  chosen  so  as  to  have  It  travel  parallel  to  the 
magnetic  field  of  the  earth  over  a  considerable  part  of  its  orbit  and  equipped 
for  regulating  its  own  charge  would  be  useful  for  the  studies. 


-12- 


References 

1.  S.  F.  Singer  and  E.  H.  Walker,  "Plasma  Compression  Effects  Produced  by 
Spava  Vehicles  in  a  Magneto  Ionic  Medium",  (Included  as  part  of  this 
report;  see  table  of  Contents). 

2.  W.  K.  H.  Panofsky  and  M.  Phillips,  Classical  Electricity  and  Magnetism. 
Addison-Wesley,  Cambridge,  Massachusetts,  (193^)* 

5.  L.  Spitzer,  Physics  of  Fully  Ionized  Bases.  Interscience  Pub  1 ishers, 

New  Ycrk,  1956. 


-13- 


Flg.  I.  Diagram  to  show  th#  Important  niiaiititles  Involved  In  the  calcjiation 
of  the  reflection  coefficient  for  an  aleetrowmgnetic  wave  incident  on 
a  layer  of  thlcknese  I.  » 


Interaction  of  West  Ford  Needles  with  the  Earth's  Magnetosphere 


Proft  S.  F.  Singer 


Abstract 

The  proposed  West  Ford  experiment  envisages  the  creation  of  a 

'belt  of  orbiting  dipoles  for  communication  purposes.  Because  of  the 

large  area-to-mass  ratio  of  the  needles  the  perturbing  effects  of 

radiation  pressure  are  of  great  importance.  It  is  argued  here  that  the 

needles  will  acquire  an  electric  charge,  and  that  perturbations  due 

to  Coulomb  drag  and  Lorentz  force  cannot  be  neglected.  Coulomb  drag 

will  lead  to  a  much  more  rapid  shrinkage  of  the  orbit  than  ordinary 

* 

drag;  their  effects  are  tabuleted.  Finally,  a  modification  of  the  West 
Ford  experinent  Is  suggested  which  could  separate  the  perturbing  effects 
this  would  lead  to  a  direct  experimental  determination  of  the  parameters 
which  enter  into  the  theory  of  the  Interaction  of  small  particles  with 
the  magnetosphere  (the  region,  from  the  ionosphere  out  to  about  10  earth 
radii,  where  the  geomagnetic  field  is  a  controlling  factor). 


1 


Introduction 

V.  E.  Morrow*  of  th«  Lincoln  Laboratory  of  MIT  has  suggested  the  use 

of  a  new  method  for  Intercontinental  communication.  It  is  based  on  the 

use  of  belts  of  orbiting  dipoles  which  act  as  a  scattering  medium  for 

microwaves.  This  technique  is  claimed  to  have  significant  advantages 

over  other  satellite  communication  techniques.  However,  it  has  been 

pointed  out  by  ethers,  particularly  observational  astronomers,  that  this 

operational  project  labeled  West  Ford  would  have  some  undesirable  side 

2  3 

effects  on  optical  astronomy  and  on  radio  astronomy.  Our  present  paper 
is  not  concerned  with  theta  effects  nor  with  communications  end  related 
topics.  Instelid  we  shall  examine  In  some  detail  certain  aspects  of  the 
interaction  of  the  West  Ford  medium  with  the  exosphere.  In  particular, 

a 

the  Interaction  of  the  electric  charge^whlch  the  needles  are  certain  to 
carry^with  the  ionized  portion  of  the  exosphere.  We  want  to  point  out 
the  possibility  of  making  unique  observations  which  can  help  establish  the 
Important  parameters  of  this  interaction.  Therefore,  as  far  as  possible, 
this  paper  will  give  numerical  predictions  of  orbit  changes,  in  addition, 
a  useful  modification  of  the  experiment  is  suggested. 

Forces  acting  on  small  particles  in  near  space 

The  unique  properties  of  the  orbiting  dipoles  arise  from  the  fact 
that  they  are  so  small.  They  have,  in  fact,  close  to  the  largest  area>to-mass 

A 

r'atio  (50  cmVgm)  of  any  artificial  satel  I  i  te^and  approach  an  area-to- 
mass  ratio  of  small  (a«100  micron)  Interplanetary  dust  particles.  Hence 
we  can  apply  to  the  West  Ford  medium  a  similar  analysis  used  for  inter¬ 
planetary  dust  in  the  vicinity  of  the  earth. ^ 

The  motion  of  a  particle  Is  determined  by  the  forces  which  act  on 


It.  within  the  region  of  the  earth's  exosphere  (taken  to  extend  out  to 
ebout  10  earth  radii)  the  gravitational  force  of  the  earth  predominates 
and  will  be  orders  of  magnitude  greater  than  other  forces  for  small 
particles.  However,  there  are  three  other  forces  which  need  to  be  taken 
into  account,  particularly  for  orbiting  particles:  1.  radiation  pressure; 

2.'  the  Lorentz  force  produced  by  the  earth's  magnetic  fieldj  and  3.  Cou¬ 
lomb  force  produced  by  the  Interaction  with  the  Ionized  exospheric  gas. 

1.  The  radiation  pressure  force  is  chiefly  due  to  the  solar  radi¬ 
ation  and,  to  a  smaller  extent,  to  the  radiation  scattered  from  and  emitted 
by  the  earth.  The  dynamic  effects  of  the  radiation  pressure  force  have 

been  Investigated  by  Jones  and  Shapiro.^  The  effects  predicted'for  the 

4 

West  Ford  medium  are  given  in  the  paper  by  Morrow  and  Made! Ian  and 

^  a 

would  indicate  that  the  medium  would  disperse  at  a  certain  rate;  also  that 
the  change  in  eccentricity  produced  by  a  resonant  phenomenon  would  be 
particularly  strong  for  polar  orbits  at  altitudes  of  3700  km. 

^  These  predictions  concerning  the  behavior  of  the  dipoles  must  be 
accepted  with  some  caution  since  they  are  based  on  a  certain  value  of 
area-to-mass  ratio,  i.e.  50  cm  /gm  and  the  implicit  assumption  that  this 
ratio  remains  constant.  Now  in  fact  there  will  be  powerful  torques  act¬ 
ing  on  these  particles  due  to  the  interaction  of  their  charge  ahd  possible  magnetic 
moment  with  the  plasma  as  well  as  with  the  earth's  magnetic  field.  This 
may  give  rise  to  a  variable  area-to-mass  ratio  which  may  be  related  to 
position  within  the  orbit.  Thus  it  Is  conceivable  that  when  viewed  from 
the  sun  this  will  lead  to  other  resonances,  and  It  is  not  Immediately 
obvious  whether  they  will  reinforce^or  cancel  and  overcome  the  resonances 
calculated  for  a  spherical  particle. 

2.  The  Lorentz  force  acts  perpendicular  to  the  particle's  velocity 
vector.  It  will  be  most  pronounced  in  the  polar  regions  where  the 
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magnetic  field  is  strong  and  always  perpendicular  to  the  particle's 
velocity  vector.  Only  one  general  statement  may  be  made  about  the 
effects  of  the  Lorentz  force:  it  cannot  change  the  total  energy  of  the 
particle.  In  other  words,  the  semi-major  axis  of  the  orbit  will  remain 
constant  but  the  other  parameters  of  the  orbit  mey  vary. 


3.  The  Coulomb  (electrostatic)  drag  force  has  been  shown  to  be  one 
of  the  most  important  forces  for  particles  moving  in  interplanetary 
space. ^  This  electrostatic  drag  force  acts  generally  between  any  charged 
body  moving  through  a  plasma  and  has  therefore  been  also  cons i dered 'for 
the  case  of  earth  satellites  by  various  authors  who  however  arrj.ve  at 
varying  conclusions.®’^’ ^ Our  view  agrees  most  closely  with 

g 

that  of  Lehnert  who  first  showed  that  the  Coulomb  drag  would  not  be 

particularly  effective  for  earth  satellites  of  nominal  dimensions. 

The  importance  of  the  Coulomb  force  then  hinges  on  two  points: 

1.  the  small  size  of  the  body,  and;  2.  the  question  of  its  charge. 

In  general  we  may  say  that  if  the  body  is  very  much  smaller  than  the 

Debye  screening  distance,  then  the  question  of  the  Coulomb  drag  should 

be  investigated  in  detail.  The  dynamics  of  the  charging  situation  was 

investigated  by  the  writer,  who  also  discussed  the  competition  between 

photoelectric  emission  and  accretion  of  charge  by  ambipolar  diffusion 

from  the  surrounding  plasma.  It  was  concluded  that  in  Interplanetary  space 

7  14. 

dust  particles  would  most  likely  be  positively  charged'.  Whipple  'has 
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pointed  out  explicitly  that  the  charge  becomes  less  positive  as  one 
approaches  the  earth,  and  becomes  negative  in  the  ionosphere. 

Since  the  Coulomb  drag  produces  a  dissipation  effect  akin  to  fric¬ 
tion,  It  will  decrease  the  semi-major  axis  of  the  orbit  continuously  and 
eventually  cause  the  particle  to  enter  the  dense  atmosphere  of  the  earth. 
We  win  therefore  be  particularly  Interested  in  Investigating  the  effect 
of  Coulomb  drag  on  the  particle's  lifetime. 

Analysis  of  the  Coulomb  drag  force.  For  the  purposes  o,f  our 
discussion  we  will  adopt  a  treatment  introduced  by  Chandrasekhar  and 
Spitzer  in  which  the  orbiting  particle  will  be  referred  to  as  the  "test 
particle^'  and  the  ions  and  electrons  of  the  exosphere  as  the  "field  par- 
‘tlcles".  The  drag  deceleration  is  given  by*^ 


where  the  subscript  I  always  refers  to  the  field  particles; 

X  X  (1 .5)  ’  w/Cj 

where  Cj  x  (3kTj/mj)^  is  the  RMS  velocity 

w  s  test  particle  velocity  ( km/sec) 

G(x)  X 

^(x)  X  2rr^  exp(-y^)  dy 
o 

4  2  2  2 

A,  X  Stt  e  n,  Z,  Z  Jin  A/m  • 

,  D  .  11 

A  »  3(2Z  Z^  e^)’’  (k^T^/Trn,)i 

where  n  end  Z  denote  particle  concentration  and  charge,  respectively^for 
* 

the  exosphere  Z^  :  I . 

1 

>  With  an  exospheric  temperature  Tj  of  1500^  K;  we  obtain:  Cj  x  6.1  km/sec 


for  protons;  Cj  ■  260  km/sec  for  electrons.  Under  these  conditions 
'>’20:  G()()  r«0.2  for  the  field  protons  and  loss  than  0.015  for 
flaid  electrons.  Since  the  particle  mass  m  »  m^ ,  Eq.  (1)  reduces  to 

aj.  1 .23  X  10“^3  2^  ^2) 

In  the  derivation  of  Eq.  (2)  some  approximations  have  boon  made; 
e.g.  the  particle's  velocity  w  was  assumed  constant  and  appropriate  to 
a  circular  orbit  at  an  altitude  of  9^00  km.  No  greater  refinement  appears 
necessary  In  view  of  the  gross  approximations  which  must  be  made  In  the 
calculation  of  the  charge  Z. 

Calculation  of  electric  charge 

The  electric  charge  of  the  particle  Is  given  by  Its  electric 
capacity  multiplied  by  the  potential  of  Its  surface.  For  a  conducting 
sphere  of  radius  a  cm 

* 

2  z  6.95  X  10^  aV  (3) 

where  V  Is  measured  In  volts.  For  a  body  of  different  shape  the  capacity 
must  be  calculated;  In  the  absence  of  detailed  data  we  will  adopt  a  sys¬ 
tem  of  two  spheres  each  with  “effective"  radius  of  1.5  x  10  ^  cm.  Thus 
for  the  West  Ford  needles  we  take 

Z  'v  2  X  10^  V  (4) 

The  calculation  of  the  potential  V  also  presents  special  problems 

t 

since  It  depends  on  the  position  of  the  particle,  and  therefore  In  turn 

! 

on  the  characteristics  of  the  particle's  orbit.  It  also  depends  on  the 
work  function W'and  photoelectric  yield  i  of  the  particle  material,  its  shape 
and  its  attitude  towards  the  sun.  As  explained  earlier,  we  calculate  V 
in  an  analysis  involving  competition  between  the  photoelectric  effect  and 
ambipolar  diffusion.^  The  ambledt  electron  density  above  1 .5  R.  (^^cth 

r.  C 
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radtl)  is  given  by  Hailiweil  as 

n,  ~  10^  r‘3  (5) 

where  r  is  measured  in  terms  of  earth  radii. 

The  photo  emission  Is  determined  by  the  ultraviolet  spectrum  of  the 
sun;  I  ,e.  ,by  Ny. the  photon  fjux  per  cn^/sec  with  energy  h>r>.  eV  h-  W. 

f-S  f  19 

V  ;  where  fr  10  and  ^  s  4.  On  the 
basis  of  HInteregger's  measurements^^  •  better  set  of  values  is 

s  2.2  X  lo'^  tcm^  se3’’  (6) 

The  equilibrlufii  potential  is  given  by^^ 

(7) 

when  the  particle  has  a  large  positive  potential  and  the  ion  current  can 
be  neglected.  In  the  earth's  shadow,  where  photo  emission  is  negligible, 

the  equilibrium  potential  is  given  by  * 

\ 

V.  =  -  (kT,/2e)  in  (mjT^/mJ,)  (8) 

in  view  of  the  uhcertaintles  In  7  ,  W,  attitude  and  orbit  orienta¬ 
tion  we  have  decided  to  make  some  drastic  approximations  and  represent 
the  potential  by 

V  =  7.5(r  -  2)  (9) 

As  can  be  seen,  V  varies  from  -7*5  volt  in  the  F-Iayer  to  '^60  volt  at 
10  ft£.  Note  that  V  =  0  at  2  Rg;  the  implication  will  be  discussed  below. 
It  should  perhaps  be  pointed  out  explicitly  that  another  powerful 

1 

charging  process  exists,  namely  the  high  flux  of  fast  electrons  in  the 
outer  radiation  belt;  it  is  however  quite  localized  and  confined  mainly 
to  the  vicinity  of  the  lines  of  force  which  extend  to  3-5  R^  in  the  plene 
of  the  magnetic  equator,  its  effect  is  to  produce  (or  enhance)  the 
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negative  charge  (by  accretion  and  secondary  electron  emission);  this  effect 
must  be  considered  for  orbits  which  intersect  the  outer  radiation  belt. 

By  combining  Cqs.  (2),  (4),  ($}  and  (9)  we  obtain  for  the  Coulomb 
drag  deceleration 

a- a-^- 3  X  10“^  (r-2)^  r‘Vm 
w 

4.-4 

With  the  needle  mass  given  as  10  gm 

ag*>w  3  X  10  ^  (r-2)^  r  ^  |^/soc^3  (10) 


This  acceleration  may  be  compared  with  the  accelerations  due  to 
(I)  Radiation  pressure 


4  X  lO^erQ/cmVsec 


RP  Vs  X  lo'®  cm/sec 


a  2.5  X  10*^  cm/sec" 


(11) 


(II)  Gravity 

ig  a  980  r*^  (12) 

(III)  Magnetic 

5  -4 

With  the  circular  velocity  w 8  x  10  r  ‘  cm/sec,  and  with  the  magnetic 
field  B'''0.5  r'^  gauss,  we  obtain 


a  a  Ze  w  B/m 
H 

a  10  ^  (r-2)  r  cra/sec^ 


(13) 


Effects  of  Coulomb  Drag 

The  orbit  changes  and  lifetime  of  a  satelTIte  in  the  presence  of  a 
dissipative  drag  force  can  be  treated  most  simply  by  energy  considera¬ 
tions.'*^ 

The  energy  loss  per  orbit  in  a  circular  orbit  is  ^E  c  2'\t  r  Fq 
and  this  will  reflect  Itself  in  a  lowering  of  the  orbit  by  an  amount 
Ar  in  a  time  which  is  simply  the  orbital  period 
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T  »  Zrr  (14) 

Since  the  total  (kinetics  potential)  energy  Is  given  by  E  a  ^  GHgtn/r, 

we  obtain  by  differentiation 

2 


Ar  a  (-aE) 


2  r 
CM  m 


(15) 


Wa  can  therefore  obtain  an  expression  for  the  rata  of  shrinkage  of  a 
circular  orbit  as  follows 


s  2tc  r  Fp 


a  i  T 
^  m 


2  r* 
GH  m 


(GM) 


1/2 


3/2 


2Tt  r 


^  (GM)  " 


(16) 


which  Is  (1/tc)  times  the  orbital  period  times  drag  deceleration. 

Ue  can  also  express  (^r/ &t)  as  a  function  of  r  using  our  Eq.  (10) 

372 

for  Coulomb  drag  and  using  T  s  5080  r  for  orbital  period: 

'  (Af/At)^s  -w5  X  10  ^  (r-2)^  r"^^^  |cm/secj  (17) 

This  expression  may  now  be  compared  to  the  neutral  drag  acceleration 
evaluated  from 

F^m  ■  i  Cp  A/  w®/m  ■  i  Cp  A(nj^  +  n^  W*/m  ( I8) 

Using  ■  (7«9  ^  IO*)*/r;  m  ■  lo"*  gm;  A  ■  5  lO”®  cit^;  C  2;  we  obtain 

F^m  -  5  X  lO"*  X  1.66  X  lo'**  (n^^  +  16  np)(6.25  x  10^^)  lOVr 

-  5.2  X  10"“  (n^^  +  16  nQ)/r  [cm/sei:;*^  (18') 

In  the  case  of  a  particle  belt  it  is  more  relevant  to  calculate  (A-r/At)/ 
rat^ar  than  the  change  In  orb  I  tel  period. 

(ir/dt)^»  l/ir  (5080  r  ^*)  5.2  x  lo”*^^  "(p/’’ 

«  8.4  X  10“®  (Oj^  ♦  16  np)  r  fcm/se<Q 


(19) 
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It  Is  of  considerable  interest  to  establish  the  relative  importance  of 
the  two  drag  forces.  In  order  to  do  this  we  must  adopt  a  model  for  the 
structure  of  the  magnetosphere. 

Observational  information  Is  limited  to  rather  low  altitudes.  The 
distribution  of  mass  density  below  300  km  Is  known  from  analysis  of  satel¬ 
lite  drag^*.  The  Integrated  thickness  of  the  neutral  hydrogen  cloud  around 
the  earth  Is  known  from  analysis  of  the  profile  of  solar  hydrogen  Lyrnan-^®. 
Electron  densities  up  to  15^  are  known  from  rocket  measurements®®.  The 
electron  density  to  several  earth  radii  is  Inferred  (with  some  uncertainty) 

I 

from  whistler  observations*^'®*. 

On  the  other  hand  the  relat Ive  distributions  of  the  neutral  con¬ 
stituents  can  be  deduced  from  a  theory  of  the  exosphere*®.  The  neutral 
components  describe  ballistic  orbits  without  collisions  above  the  base 
of  the  exosphere  ($30  km)  and  the  distribution  of  concentration  with 
altitude  Is  calculable  In  a  fairly  straightforward  manner.  The  ioni,zed 
components  (having  much  larger  cross  section  because  of  Coulomb  interaction) 
are  distributed  according  to  the  barometric  formula®^.  In  each  case  the 
slope  (or  scale  height)  Is  determined  by  the  temperature  at  the  base  of  the 
exosphere.  Using  the  data  referred  to  above  for  normalization,  it  is 
possible  to  construct  a  model  of  the  exosphere*®.  This  model  is  shown 
In  Table  I. 
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TABLE  I 

Concentration  of  major  constituents  In  the  terrestrial  exosphere  (cm”^) 


r/R 

0 

0* 

H 

l.lOO 

1.2  X  10^ 

5  X  10® 

10  X  10® 

1.200 

2.7  X  10* 

2.5  x.lO* 

6.6  X  10® 

9x10® 

1. 300 

1 .7  X  10^ 

1.7  X  10® 

4.4  X  10® 

6.6  X  10® 

I.LOO 

1 

1.3  X  10^ 

3.1  X  10* 

5  X  10® 

1.500 

— 

II 

2.3  X  10® 

3.7  X  10® 

1.75 

<.I 

1.25  X  10® 

2.1  X  10® 

2.00 

8x10* 

1.3  X  10® 

3.00 

2  X  10* 

‘  3*5  X  10® 

4,00 

82 

1.5  X  10* 

5.00 

43 

85 

6.00 

i 

25 

51 

7.00 

16 

55,. 

8.00 

12 

23 

9.00 

9.1 

17 

10 

6.5 

11 

n 


TABLE  2 


Drag  Deceleration  and  Orbit  Shrinkage 

(i  r/At) 

a 

c 

r/R 

( cm/sec®) 

( cm/sec) 

( cm/sec®) 

(cm/sec) 

I'l 

9  X 

10'» 

1.7  X  10 

1.8  X  10“» 

3.5  X  10"® 

1.2 

3.6  X 

1 

O 

7.6  X  lo"* 

1.1  X  10“® 

2.4  X  10“® 

1.3 

1.3  X 

10"« 

3.6  X  lo"* 

6.9  X  10“® 

1.7  X  10*® 

1.1* 

3.7  X 

10*^ 

1  X  lo"® 

4.0  X  I0“® 

1.1  X  lo"® 

1.5 

2.1  X 

10“7 

6.1  X  10““* 

2.2  X  10“® 

6.8  X  10“* 

1.75 

9.8  X 

lo'* 

3.7  X  10“* 

3.6  X  10“^ 

1.4  X  10"* 

2.0 

5.5  X 

10"* 

2.5  X  io"* 

(0) 

(0) 

3.0 

9.5  X 

• 

1 

o 

7.9  X  10“* 

1.1  X  I0“® 

9.6  X  10"* 

4.0 

3  X 

10’» 

3.9  X  10"® 

1.9  X  10"® 

2.5  X  10"* 

5.0 

1.4  X 

10’» 

2.4  X  10"® 

2.1  X  10“® 

4.0  X  10"® 

Both  nautr«]  and  Coulomb  drag  dacalaratlons  ara  shown  In  Table  2, 
as  wall  as  the  resulting  values  of  (Ar/dt}t  It  must  be  cautioned  however 
that  the  Coulomb  drag  depends  on  many  detailed  factors,  (some  of  which 
were  discussed  earlier),  involving  the  particular  properties  of  the 
particle,  of  the  exosphere,  and  of  the  sun;  as  well  as  the  particulars  of 
the  orbit. 

The  zero  value  of  Coulomb  drag  at  2Rg  should  not  be  taken  literally 
it  merely  Indicates  that  somewhere  in  this  altitude  range  there  may  be 
a  minimum.  As  discussed  elsewhere,  small  fluctuations  of  charge  will 
produce  a  drag  (since  it  depends  on  the  square  of  Z  and  not  on  the  sign); 
In  particular,  the  day-night  effect  sets  a  lower  limit  to  the  Coulomb 
drag^.  This  m/nimum  may  give  rise  to  a  high  concentration  of  inter¬ 
planetary  smoke  (particles^  1  micron)  near  2  earth  radii  if  they  exist. 
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Discussion 

It  would  be  of  Interest  to  calculate  the  lifetimes  of  the  West  Ford 
particles.  Since  the  details  of  the  orbit  have  not  been  revealed  and,  In 
particular,  since  the  altitude  and  orbit  inclination  Is  not  known. 

It  does  not  seem  worthwhile  to  attempt  any  refined  calculations.  There 
are  in  any  case  uncertainties  about  the  electric  charge  of  the  particle, 
about  Its  attitude  with  respect  to  magnetic  field  and  with  respect  to  the 
sun,  and  also  about  the  actual  values  of  air  density. 

With  Initial  orbit  details  available,  the  subsequent  history  of  the 
West  Ford  particles  can  only  be  obtained  by  numerical  methods.  It  Is 
clear,  by  referring  to  Table  2,  when  the  neutral  drag  may  be  neglected. 

It  is  clear  also,  by  referring  to  £qns.  (10),  (11),  and  (1?),  that 
probably  none  of  the  forces  can  be  neglected,  except  for  very  special 
orbits. 

In  a  numerical  computation,  particular  attention  must  be  paid  to 
the  day^nlght  effect,  not  only  on  radiation  pressure,  but  also  on  the 
particle  charge.  In  fact,  because  of  fluctuations  in  radiation  belt 
intensity,  solar  flares,  and  similar  factors,  the  charge  will  fluctuate 
unpredictably.  This  effect  will  contribute  to  a  spreading  of  the  belt 
which  is  in  addition  to  the  dispersion  calculated  by  the  Lincoln  Laboratory^. 

^  It  is  interesting  to  note  the  importance  of  the  radiation  pressure 
in  relation  to  the  Coulomb  drag.  This  is  due  to  the  geometric  shape  of 
the  needles. 

it  should  be  noted  that  for  a  spherical  interplanetary  dust  particle 
the  radiation  pressure  acceleration  is  x  lO’^  Tra^/m  and  Is^relatively^ 
100  .  1000  times  smaller  when  compared  to  the  Coulomb  drag^'*^. 


I 
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NoverthelesSj  the  shrinkage  of  an  Initially  circular  orbit  produced 
by  energy  loss  due  to  drag  Is  a  gross  effect  which  should  be  observable^ 
and  which  gives  some  Information  on  air  density^  as  well  as  about  the 
Importance  of  Coulomb  drag. 
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A  Proposed  Experiment 

However,  a  much  more  refined  experiment  is  possible  which  illuminates 
the  interaction  of  small  particles  with  the  magnetosphere.  It  is  based  on 
a  separation  of  the  perturbing  forces. 

Consider,  for  example,  an  initially  circular  orbit  at*^3200  km 
altitude  (r  >  I.5  R^),  nearly  polar  and  perpendicular  to  the  earth-sun 
line.  In  this  special  orbit  the  radiation  pressure  does  not  produce  any 
appreciable  perturbation  (at  least  Initially),  and  the  Lorentz  force 
perturbation  can  be  determined  quite  simply,  although  it  is  again 
negligibly  small.  Since  the  particle's  velocity  vector  Is  moreor  less 
in  the  plane  of  the  magnetic  meridian,  the  chief  perturbation  will  be 
orthogonal  and  cause  a  rotation  of  the  orbital  plane.  The  sense  of  the 
rotation  establishes  the  si  on  of  the  electric  charge  on  the  particle. 

Using  the  analysis  leading  up  to  Eq.  (I3)  we  estimate  the  rate  of  rotation 

of  the  orbital  plane  as  about  rad/sec. 

i 

Next  consider  a  modification  of  the  photoelectric  properties  of  the 

•  a 

copper  needles.  Specifically,  one-third  could  be  left  with  their  present 
copper  surface,  one-third  would  be  coated  with  a  material  of  high  work 
function  W  and  low  photoelectric  yield  (e.g.  platinum)^  and  the  remainder 
with  a  material  of  opposite  characteristics  (possibly  a  cesium  alloy). 

For  our  particular  orbit,  the  Cu  particles  would  be  negative,  the  Pt 
particles  even  more  negative,  while  the  Cs  particles  would  have  a 
nearly  zero  potential,  or  perhaps  even  positive. 

Because  of  Coulomb  drag  the  three  belts  will  shrink  at  a  different 
rate.  uMhg  the  analysis  leading  up  to  Eq.  (17)  we  estimate  a  differential 
rata  of  shrinkage  of  orbit  radius  amounting  to  nearly  I  km  par  day. 


Thus  the  separation  of  this  West  Ford  belt  Into  three  layers  might 
be  resolvable  with  the  available  radar  methods,  provided  the  dispersion 
of  particles  at  injection  can  be  held  very  small. 
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Further  on  the  Interaction  of  West  Ford  Needles 
vlth  the  Earth's  Magnetosphere. 

S.  F.  Singer,  Riyslcs  Depeurtment 

university  of  Maryland 
College  Park,  Maryland* 

OSie  purpose  of  the  present  paper  Is  to  amend  and  expand  my  earlier 

discussion  of  the  "Coulomh  drag"^.  Originally  the  Interest  lay  In 

calculating  the  drag  force  a(^lng  on  a  charged  spherical  duat^..p§;|;^lc3.e 

2 

moving  through  an  Ionized  gas  ;  the  treatment  can  be  extended  to  the 

3 

case  of  the  orbiting  dipole  needles,  used  In  Project  West  Ford  .  Now 

that  much  of  the  xmderlylng  theory  has  been  published  and  the  orbit  has 
4 

been  detailed  ,  it  becomes  possible  to  discuss  the  Ixifluence  of  Coulomb 

drag  on  the  lifetime  of  the  needles. 

4 

Details  of  West  Ford  Experiment 

Ihe  copper  needles  have  a  length  J?  of  1.77  cm,  <uid  a  diameter  d 

.  -4  " 

of  2.0  X  ID  ^  cm.  The  mass  H  is  10  gm.  They  are  supposed  to  rotate 

in  propeller>llke  fashion  at  about  2  rev.  per  second;  their  (projected 

area)/mass  ratio  A/m  will  vary 


*  Present  address:  Jet  Propulsion  laboratory,  California  Institute  of 
Technology,  Pasadena,  California 
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betveen  a  maximum  value  of  $0  cm^/gm  and  a  minimum 
value  of  3^  cm^/gm. 

Obe  Initial  orbit  Is  supposed  to  be  a  nearly 
circular  polar  orbit  (l.e.  with  inclination  *6  >  $0° 
and  eccentricity  e  /w  0 }  at  an  altitude  of  about 
3800  km  above  aealevel^or  with  semlmajor  axis  CZ.  ^ 
1.6  R  .(earth  radii)  from  the  earth's  center,  ^e 
orbital  plane  is  supposed  to  make  a  small  angle 
itiltlally  with  the  sun-earth  line. 

Radiation  Pressure  Resonance  Effect 

For  special  orbits,  and  for  a  certain  band 
of  orbited  parameters  near  these  orbits,  there 
e::lsts  an. effect  produced  by  solar  radiation  pressure 

vhlch  can  limit  the  lifetime  of  a  satellite.  For  a 
resonant  orbit  the  eccentricity  e  will  Increase,  with 
no  appreciable  change  in  a;  essentially  the  orbit 
Is  displaced  within  its  plane  until  the  perigee 
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reacliea  the  deneer  Rtmuspltnrr  lutd  t.ho  orbit  ulirinkfi 
and  decays.  Tills  is  the  picture  as  presented  by 
Shapiro  and  Jonea^  Obey  shov  that  the  lifetime  of 
the  West  Ford  needles  vlU  be  about  7  years  for  the 
Initial  orbit  chosen.  !Qie  range  of  the  resonance 
Is  rather  narrow;  for  kjn  -  $0  cm  /gm,  the  altitude 
range  is  ±  ISO  km,  the  range  of  Inclination  Is  ±  1°. 
For  a/m  ■  33 >  the  resonance  region  Is  considerably 
narrower  but  not  quoted;  the  lifetime  Is  increased 
i  In  Inverse  proportion  to  a/n,  l.e.  to  about  10 
years. 

Without  detailed  numerical  orbit  calculations 
It  seems  difficult  to  establish,  e.g.,  the  lifetimes 
near  the  edge  of  the  resonance  region^ or  its  width. 
However  we  can  quite  easily-  establish  the  dependence 
of  CL  ^or  the  resonant  orbit  on  1.  Qie  rate  of 
change  of  the  argument  of  a  nearly  circular  orbit 
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perlgee  1b  given  by  (see  e.g.  Ref  $) 

ci)  '0  (1) 

Oxe  regression  of  the  node  of  the  orbit  Is  given  by 

jj_  A,  tC  (^/r)  (2) 

and  oust  be  added  to  ^^the  rate  of  rotation  of 
the  sun  with  respect  to  the  Initial  orbital  plane. 

Of  course,  ^^1.  (All  angular  rates  are  given  In 
degrees  per  day).  The  resonance  condition  demands 
-  cL  /  (3) 

For  1  near  90^,  and  for  a'^  1.6  R  It  foUows  that 
-C  j»4(r  decree 

Aaavmlng  that  the  needles  are  launched  Into  a 
perfect  resonant  orbit,  ve  must  Investigate  If  there 

k 

:.exlst  aj^  perturbations  which  can  spoil  the  resonance. 
For  example,  a  drag  force  which  decreases  CL  by  about 

$0  km  par  year  would  be  sufficient.  In  a  matter  of  .only 

3  years  to  pull  the  orbit 

/ 


i 
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out  of  resonance,  provided  the  inclination  is 
cbaxtged  at  the  same  time  by  a  lesser  amount  than 
given  by  Sq.  (4). 

Ifeutral  drag  is  clearly  too  small,  leading  to 
a  rate  of  change  •mci,  ^  .  0.1  km  per 

year^.  Ohe  Coulcnb  force  appears  to  be  the  only 
one  of  iagKMrtanee .  It  vlU  be  examined  quantitatively 
later;  let  it  be  denoted  by  F.  Ve  must  note  that 
no  matter  how  its  magnitude  may  fluctuate,  it  will 
always  be  directed  opposite  to  the  relative  velocity 
vector  between  the  needle  and  the  plasma.  We  consider 
the  latter  to  be  corotating  with  the  earth's  angular 
velocity^  at  least  as  far  up  as  aurorsd  latitudes, 
say  A  ^  6$^.  Hence  the  needle,  assumed  to  be  in 
a  polar  orbit,  will  experience  a  cross  conqponent  of 
force  F,^n  0  A 

where  ■ioh  e  =  '■/x;!  “  «»  ' 
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For 


JT  Jr(.4»  p'  J  ^  "hxv^O  &  ~  O.IZ 


We  now 


the  ehrlakege  of  the  orbit. 


given  by 

(Ai-Zcwt)  -  -(TA)(;F/M),_(zr/v)(F/M)  (5) 


where  T 


iB  the  orbital  period  i^/k 


During  an  Interval  the  Increment  In  croaa 
velocity  la  ^ 

AVi-  At/r  /([P/m)  -sivv  ^  /Cci  X  <it. 
c 


(6) 


Oiercfore  the  change  in  inclination  vlU  be  given 
by  Ac  •  (aVj^/V/)  (7) 

We  can  now  calculate  the  reaultant 

/Uv  k»vv  yet»- .  (®) 

Note  that  thla  ratio  la  Independent  of  F>  For  r 

1.6  R 


Ao./^  A.  4'.  8  X 


(o 


Ion  per  degree.  (9) 
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oilt  result  is  about  50  times  greater  than  the 
slope  of  the  "llae  of  resoaaaoe",  given  by  Sq.  (h). 

Ve  conclude  therefore  that  a  drag  force  yleldli^ 
( Cu  [  >  20  km  per.  year  will  Indeed  pull  the  orbit 
out  of  resonance.  Ve  therefore  arrive  at  the 

paradoxical  result  that  an  increased  drag  could 
lead  to  a  much  longer  lifetime.  However,  if  the 
drag  is  very  much  greater,  say  1  ^  |  >’ 

500  km  per  y«'ar,  then  the  lifetime  could  be  even 
shorter  than  the  value  of  7  years  deduced  from 
resonance  considerations. 
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Calculatlon  of  Coxtlomb  Drag  for  a  Needle 

Tb«  CAloulatloa  of  tbo  dreg  force  coaeiste  of 
three  steps:  (l)  the  electrical  capacity  C  ; 

(ii)  the  electrical  potential  ,  and  finally 
(ill)  the  drag  force  f"  . 

(l)  Capacity.  In  spite  of  the  admonition^  that 
"for  a  body  of  different  shape  jtTtm  a  spher^  the 
capacity  must  be  calculated"^  I  omitted  to  do  this 
and  adopted  an  incorrect  value  which  turns  but  to 
be  appropriate  to  a  sphere  of  diameter  i[z 
cm  but  not  to  a  needle.  If  we  approximate  the 
dipole  by  a  prolate  spheroid  (as  suggested  to  me 
by  Professor  leverett  Davis),  its  capacity  is 
given  (in  MKS  units) 

C  •  4tc  e  (qJ-  b’-/*’  [tanfi  ‘  ((-  (lO) 

r* 

while  a  sphere  of  radius  a  has  a  capacity 
=  4'^d:  CL  ^  !•  112.8  KIQ  CL 


(11) 
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A  lavet  limit  to  C,  the  actual  capacity  of  the 
needle.  Is  obtained  by  taking  ct* 
and  b  * 

Since  h  «  CL  ,  ve  can  approximate 

Cj,,  ATcecL^i  (12) 

'''  4tc6  cl  C^-A-sJ  ' 

vhere  (  J  -  -  I  -  i 


An  upper  limit  to  C  la  obtained  from  the 
capacity  of  the  apheroid  which  clrcumacrlbea  the 
needle,  assumed  to  be  cylindrical;  =  \fz  C^r 

I  therefore  adopt  a  mean  value 

^  ^  “O.lJt  'v  2.0  u.  (^3) 

It  la  interesting  to  note  that  the  sphere  \dilch 
circumscribes  the  cylinder  has  a  capacity^  $.5 

t 

times  that  of  the  needle* 

For  the  case  of  the  West  Ford  needle, 

C-o.t«  C±2o7.j  (14) 


/>n/ 
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Zbe  quoted  iaeccurecy  eteoe  partly  from  our  Ineufficlent 

knowledge  of  the  detailed  shape  of  the  needles 

(straight  or  bent?  ends  flattened?  smoothed?).  It 

should  be  noted  that  O  Is  about  ^8  times  greater 

than  my  previously  adopted  value  in  Ref.  1;  hence 

the  needle  charge  *Z.  should  be  increased  by  the 

2. 

same  factor,  and  the  drag  ly  (^)  . 


-11- 


(11 )  Electrical  Potential.  As  stated  before,  the  potential  V, 
and  therefore  the  number  of  electronic  charges  Z,  are  the  most  uncertain 
quantities  In  the  calculation  of  the  drag.  Depending  on  the  position 
of  the  body  In  space,  three  processes  are  of  Importance:  (a)  accretion 
of  Ions  and  electrons  from  the  magnetospherlc  plasma;  (b)  ejection  of 
photoelectrons  by  solar  ultraviolet  radiation;  (c)  absorption  of  high 
velocity  electrons  (and  protons)  trapped  In  the  earth's  magnetic  field. 
Charging  Rrocesses:  Before  attempting  any  numerical  estimates  for  the 
West  Ford  orbit.  It  Is  of  Interest  to  make  some  general  remarks  on  the 
calciilatlon  of  charging  rates  for  the  three  processes.  Both  (a)  aiid 
(b)  depend  on  the  body's  potential  V,  vhlle  (c)  Is  reasonably  Independent 
of  V;  as  a  result  (b),  and  more  particularly  (a),  are  "compensating" 
processes  vhich  adjust  their  rates  In  such  a  way  as  to  minimize  V; 

(process  (a)  can  even  change  its  sign).  Hie  basic  reason  for  this  is 
that  (c)  Involves  energies  much  greater  than  V,  while  (a)  and  (b)  do 
not. 

Next,  It  should  be  noted  that  for  process  (c)  the  omni  directionally 
averaged  area  Is  relevant  (assuming  that  the  trapped  particles  are 
Isotropically  Incident).  For  (b)  the  unldlrectionally  (l.  e.  along  the 
sun  line)  averaged  projected  area  matters,  while  for  process  (a)  the 
aiea  depends  on  V  In  a  complicated  way,  as  well  as  on  the  shape  of  the 
body,  orientation  of  the  magnetic  field,  and>  temperatiire  of  the  plasma. 
(Oils  problem  Is  being  Investigated  in  a  fundamental  way  by  Mr.  E.  H.  Walker 
at  the  Uhiversity  of  Maryland  as  part  of  his  Ph.D.  dissertation.) 

Superiaqposed  on  these  fundamental  problems  are  more  detailed  ones 
arising  from  uncertainties  of  the  jhyslcca  parameters: 
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(a)  Of  crucial  loqportance  is  the  velocity  distribution  of  the 
plasma  electrons,  particularly  the  existence  of  a  high  energy  tall. 

The  Ions  make  good  thermal  contact  vlth  the  neutral  atoms  belov  the 
exosphere  base;  hence  they  are  reasonably  certain  to  have  a  Maxwellian 
distribution  appropriate  to  a  tesqperature  of  about  1$00°K,  l.e.  kinetic 
energy  <~0.195ev.  The  electrons  however  are  ejected  with  energies 
higher  by  1  or  2  orders  of  magnitude  and  must  thermallze  through  collisions. 
It  Is  quite  likely  that  a  goodly  fraction  of  '^e  high  energy  tall  escapes 
upward  Into  the  magnetosphere  (particularly  at  higher  latitudes  where 
the  lines  of  force  are  more  nearly  vertlccd)  for  two  reasons:  the 
path  length  Is  shorter;  and  the  Coulomb  scatterlx^g  cross  section  decreases  , 

V 

as  (l/energy)^. 

Consider  the  situation  In  the  e$urth's  shadow  and  assume  only  process 
(a)  effective.  The  equilibrium  potentlsd  would  be  given  by 

V  -  -  (kT^/2e)  la  (mi  Xg/  »e^i)  *'  (15) 

If  we  adopt  for  an  electron  "temperature"  appropriate  to  their 
non-Msxwellian  distribution, 

Te  -  (2/3k)  (Fa^)  “  ^^i  (l6) 

then  V-0.75f  volts  for  an  O’*"  plasma,  and  ^25^1  less  for  sn  ff*"  plasma. 

With  f  between  1  and  perhaps  10,  we  amy  adopt  a  pzwlslonal  value  for 
V  of  -  3*6  volts. 

At  an  altitude  of  3^00  km  we  adopt  a  density  of  3  x  10^  H'*'  Ions 

per  cm3.  Under  these  conditions  the  accretion  rate  of  positive  ions 
8 

is  '<-'10  per  eec,  very  roughly. 


(b)  For  the  flvuc  of  solar  photons  vith  energy  greater  than  hi/Ve 
had  adopted^ 

Ny,«  f  v"^  -  2.2  X  lCp-3  /"om2  sec/"^'  (1?) 

However  the  expression  Eq.  (7)  for  the  equlUbrliua  potentled  qaoted 

earlier^  applies  only  >dien  the  body  has  a  large  positive  potential,  and 

Is  therefore  Inapparoprlate  for  the  West  Ford  orbit.  Major  uncertainties 

In  process  (b)  Involve  the  photoelectric  q\iantum  yield  /  and  the  work 

function  V.  With  some  crude  estlxnates,  0.1,  W^IO  volts  we 

derive  a  photoelectric  charging  flux  of  (3  x  10^®  ca"^  sec"^)j((5  x  10“3  cm^) 
a 

<»>’  1.5  X  10^  per  see,  cooparable  to  the  flux  of  positive  Ions. 

(e)  Since  maxlnm  high  energy  proton  fluxes  In  the  Inner  radiation 
belt  only  amount  to  ^10^  per  cm^'Sec,  we  need  to  consider  only  the 
effects  of  high  energy  trapped  electrons.  Mainly  because  of  Instrumental 
considerations  we  can  divide  them  into  three  energy  categories:  (l) 

'  >  100  kev  with  a  peak  flux  of  .-^10^  per  cm^  -  sec,  located  between 
3  -  4  R  In  the  equatorial  plane,  hence  at  magnetic  latitude  50'^  near 
the  earth  (Explorer  VI  data  of  Winckler  et  al;  Fan,  Meyer,  Simpson; 

Farley  and  Rosen;  Lunik  data  of  Vemov  et  al)  (il)  "Aurored"  electrons 
of  rather  lower  energies  (10  -  50  kev)  with  flux  of  —  10^®  per  cm^-sec 

(Rputnlk  3  data  of  iOrMsovsky;  rocket  data  of  Mcllwain) 

(ill)  A  rather  hypothetical  electron  component,  assuned  to  be 
responsible  for  the  magnetic  storm  ring  current  (rather  than  the  conven¬ 
tionally  assumed  20  kev  protons).  An  average  magnetic  storm  demands  an 
energy  density  of  ^  50  kev  per  cm3;  if  the  mean  energy  of  these 
electrons  Is  300  ev  their  flux  will  be  2  x  10"  per  cm?  -  sec.  (Singer). 
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Ibe  omnidirectionally  averaged  area  of  a  needle  la  jk)  Xd 
— ^4  X  10"^  cm^.  A  mean  thickness  Is  '~'2  x  10"^  gm  -  cm"^;  hence  a 
good  fraction,  say  ^  of  cooqponent  (l)  will  be  absorbed  (if  we  Include 
the  effects  of  scattering).  As  we  go  to  lower  energies  the  effects  of 
secondary  electron  emission  as  well  as  back  scattering  become  of  increasing 
Importance.'^  Ihus  we  estimate  a  charge  accretion  of  ^  for  coo^nent  (il), 
and  only  l/lO  for  component  (ill).  ^ 

Thus  we  would  estimate  that  for  particular  portions  of  the  orbit, 

O 

constituting  perhaps  20^  of  it,  we  have  an  electron  accretion  of' — '10 

per  second,  with  a  variability  of  e  1  ~  2  orders  of  magnitude. 

For  the  West  Ford  orbit  (polar,  at  altl'tude  38OO  km).  It  would 

appear  from  the  preceedlng  discussion  that  the  charge  Is  always  negative, 

* 

being  largest  on  the  nlghtslde  In  the  auroral  zone,  and  smallest  on  the 
day  side  at  low  latitudes. 

Adopted  Values;  For  the  purpose  of  further  calculation  we  adopt  a 
potential  V  of  •  3*6  volt,  and  estimate  a  high  probability  tbltt  It  stays 
within  a  factor  of  2,  l,e, between  -1.6  volt  to  -7.2  volt. 

It  should  be  pointed  out  specifically  however  that  one  cazmot  really 
average  the  potential  over  an  orbit,  or  over  a  period  of  time.  Since 
the  drag  depends  on  V^,  It  should  be  pointed  out  explicitly  also  that 
for  a  small  particle  the  drag  is  Independent  of  the  sign  of  the  potential; 
for  a  large  body  a  negative  potential  leads  to  a  lower  drag  'than  the  same 
value  of  positive  potential.  Averaging  could  lead  'to  a  serious  under¬ 
estimate  of  the  effects  of  possible  peaks  In  the  fluctuation  of  the 
particle's  potential.  To  gauge  -these  peaks  we  calculate  -the  time 
constant  Involved  In  the  charging  process;  with  a  nonmnAi  charge  of 
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'  CV/e  ■  4  X  10^  electronSi  this  time  turns  out  to  be»'-*—0.04  seconds. 

I  am  not  aware  of  any  processes  which  have  a  shorter  time  constant  and 
would  therefore  be  effective  In  producing  peaks;  but  their  existence 
cannot  be  ruled  out  In  advance. 

(Ill)  Bie  Drag  Force;  A  massive  body  carrying  a  charge  Ze  moves 
with  velocity  w  through  a  plasma  with  Ions  of  mass  m^,  charge  Z^,  and 
of  density  ni.  dbe  drag  force  may  be  deduced  from  expressions  given 
by  3pltzez^: 

p  a  (4'?7'  e^  In  A  )  nj^  Q  (x)  a  -  6.8  x  10“^3  ^2  q  ^^.Q) 


Table  I  shows  G  (x)  for  various  conditions;  x  >  1.5^*^  (v/Ci)  idiere 
Is  the  RMS  velocity  of  the  field  particles,  corresponding  to  a  teiqperature 
^  of  1500°K.  It  can  be  seen  that  below  2000  km,  O*  Is  about  5  -  lO^t 
as  effective  as  while  the  electrons  cure  about  ^  as  effective  as  H^. 
Above  2000  km  only  needs  to  be  considered.  Seyond  10-15  B  ve  must 
consider  Interplanetary  plasma  velocities  and  tenq^ratures. 

For  a  range  of  orbits  near  the  West  Ford  resonant  orbit  we  may 

take 

F  -  1.0  X  10"23  ^2  ^^^9) 
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•  Table  I  Selected  Values  of  G(x) 


Region: 

Ionosphere 

0* 

electrons 

1000-6000  km 

2-lOR 

300-2000  km 

V  (circular) 

0.12 

m 

0.007 

.018 

V  (escape) 

0.06 

0.003 

• 

0 

ON 

It  is  useful  also  to  Indicate  the  sensitivity  of  G  to  small  charges 
In  w. 

Tor  between  IOOO-60OO  Ism  dO/dv  -  O.O3  per  km/sec; 
Tor  0*^  between  300*2000  km  -  O.OO3  per  km/sec 
Tor  electrons  *  O.OO3  per  km/sec 
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Wlth  a  nominal  charge  Z  of  4  x  10^  electrons  and  with  n^  ■  3  x  ICr 
cm"-  we  obtain  F  -  4.8  x  lO””^  dynes. 

Before  this  value  can  be  accepted,  It  Is  necessary  to  prove  that 
,  the  size  end  shape  of  the  body  is  of  minor  importance  emd  that  the  drag 
depends  only  on  Z^.  I  have  been  able  to  show  that  the  "effective  area" 
of  a  needle  is  closely  the  same  as  that  of  a  point  charge.  The  potential 

I 

<j)  in  the  vicinity  of  very  thin  straight  wire  segment  of  length  (  is 
given  by^ 

^  -  (aZe/j?)  tanh-1  (^/r^^+rg)  (20) 

where  rj_  and  rg  are  the  distances  from  the  ends  of  the  needle  to  point  F; 

the  eq]uipotential  lines  are  clearly  ellipses. 

Ihe  potential  on  the  equatorial  plane  of  the  inscribed  ellipsoid 

can  be  found  by  mailing  use  of  the  expressions 

tanh'^x  ■  sinh'^  ■ln2u  +  il  -31 

^  5  ¥u2 

i  ■ 

^re  X  ■  u  (u2  +  1)  "i. 

m  our  ease  x  ■  1  -  i  b2/e2  >  i-  (2x10”^)  u  'vjOO;  and  tanh"^  ^ 

In  2tt  >  6.9. 

We  now  wish  to  find  the  ellipse  where  the  potential  is  about  l/lO 
of  the  surface  value.  This  figure  Is  chosen  to  be  roughly  equal  to  the 
ion  kinetic  energy  and  therefore  defines  approximately  the  line  where 
the  deflection  of  the  ion  is  90°.  Hence,  setting  tanh~^  -  6.9/10,  we 
obtain  (from  ^Babies)  0}/ri  +  ^g)  .  0,5  or  rj^  +  rg  -  2.0j(.  Ihis  ellipse 
turns  out  to  have  semimajor  axis  A  •  jL,  sasiminor  axis  3  m  Z  “(V^  and 
area 

.  ^^o  -  ■  ^  *  o>\n  Ji, 


fei) 
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IBxli  area  is  to  be  compared  to  the  case  vhere  the  needle  charge 
obtained  from  Eq.  (14) 

Qfi  O.t&t  (4116  0.  (22) 

Is  concentrated  at  a  point.  Die  distance  r,,  at  vhlch  the  potential 
Is  1  V  Is  given  by 

'/o  N/  ‘  Qn)  3  i.e. 

Ohe  area  enclosed  by  this  circle  is  TQ  £/z)  (23) 

vhlch  Is  closely  equal  to  the  value  of  Eq.  (21). 

Note  that  a  calculation  for  the  two  areas  correspozuUng  to  the  full 

potential  V  gives  the  values  j^d  for  the  needle  axid  l/lOO  x  O.STlCi. 

for  the  point  charge;  they  are  in  the  ratio  of  1:17.2.  However  since  the 

contribution  to  the  drag  comes  predomlnently,  l.e.  by  a  factor  i>v  2.0 

>» 

from  small  deflections  produced  by  distant  collisions,  ve  are  Justified 
in  stating  that  the  total  charge  only,  Independent  of  body  shape',’ 
dirtemlnes  the  Coulomb  drag  force  to  a  close  approximation. 

Oils  statement  evidently  breaks  down  where  the  body's  dimension 
becomes  larger  than  either  the  Debye  length  .A/  •  6.^  ^ 

or  the  distance  r^  ,■ 
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Effect  of  Coulomb  Drag  on  Orbit 

Lifetime;  Uie  major  effect  of  Coulomb  drag  Is  a  transfer  of  energy 
to  the  plasma  in  the  magnetosphere.  Ohls  leads  to  a  decrease  in  the 
total  energy  of  the  body  and  to  irreversible  shrinkage  of  the  orbit. 
For  a  near  circular  orbit  at  distance  r  from  the  earth's  cehtbr  the 
rate  of  shrinkage  la  given  by  Eq*  (5). 


— k 

For  the  Vest  Ford  case,  vlth  M  ■  10  gm  and  a  plasma  density  at 
l.oR  of  oX  lOoMwe  obtain  the  foUovlng  result: 

^'/4't  "  16  C*»v/s*c  i  or  Sooo  km  fcr  v^ea.r  ^  or  1+ 

Dlls  result  refers  to  our  nominal  "average"  potential  of  -3*6  volt;  for 
1.8  volt  dr/dt  12^0  km  per  year;  for  -7.2  volt  dr/dt  20000  km  per 
year,  about  35  km  per  day. 

It  can  be  seen  that  the  Coulomb  drag  provides  an  upper  limit  to 
the  lifetime  of  a  few  weeks  to  a  few  years,  rather  less  than  the  life¬ 
time  obtained  from  the  radiation  pressure  resonance  effect. 

A  fuller  interesting  (and  seemingly  paradoxical)  point  exists, 
under  the  extreme  assumption  of  a  potential  as  low  as  O.36  volts  the 
rate  of  shrinkage  is  only  50  km  per  year;  it  would  pid.1  tlie  orbit  out 
of  resonance  in  about  3  years  end  thereby  increase  the  lifetime.  But 
I  estimate  the  probability  for  such  a  small  potentlcd  as  less  than  10^ « 
Dispersion;  de  dispersion  of  the  orbiting  dipoles  Is  likely  to  be 
dcminated  by  the  Coulomb  drag.  Fluctuations  in  the  chargltig  rate  are 
produced  by  the  day-night  effect,  and  less  predictably  by  solar  flares, 
magnetic  storms  etc,  Sie  resultant  dispersion  might  amount  to  a  few 


; percent  of  the  shrinkage  rate^  i.e.  to  a  spread  of  about  1  km  per  day. 

As  discussed  In  the  earlier  paper, ^  by  coating  the  copper  needles 
In  a  suitable  vay  the  work  function  can  be  altered  and  the  photoelectric 
charging  rate  can  be  modified.  One  can  conceive  of  an  experiment  In  vhlch 
three  groups  of  needles  with  different  work  functions  are  ejected;  they 
vould  decay  at  an  average  rate  of  14  km  per  day  but  with  a  difference 
amounting  to  a  few  km  per  day. 

Rotation  of  Orbital  Plane;  For  a  nearly  polar  orbit  the  Lorentz  force 

Ze  (v  X  B)  produces  an  orthogonal  perturbation  vhlch  Is  easily  estimated 

by  a  simple  geometric  perturbation  approach.  It  leads  to  a  torque  vdiich 

precesses  the  orbit  plane  at  a  rate  of  about  5  x  lO"^^  radians  per  see, 

-2 

or  1.$  X  10  rad  per  year;  It  can  be  considered  as  negligible. 


ConclusloQs 


Hhe  Importance  of  the  Couloub  drag  on  the  orbits  and  lifetime  of 
the  West  Ford  needles  is  established  beyond  doubt.  The  upper  limit 
to  the  lifetime  ranges  from  a  few  weeks  to  a  few  years;  the  chief 
uncertainty  is  the  electric  potential  acquired  by  the  orbiting  dipoles. 
Corresponding  to  a  “minimum"  potential  of  -1.8  volts  we  calculate 
a  shrinkage  rate  of  the  orbit  of  ^  4  km  oer  day;  for  a  "maximum"  potential 
'of  -7*2  vol^s  the  rate  is  ^  55  4m  per  day. 

Clearly,  therefore,  the  observation  of  the  orbit  decay  due  to 
Coulomb  drag  gives  direct  Information  on  the  average  potential  of  the 
orbiting  particle,  and  -  less  sensitively  -  on  the  density  of  plasma 
in  the  earth's  magnetosphere. 
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As  a  satellite  presumed  to  be  neutral  moves  through  the  upper  atmosphere, 
it  sweeps  out  the  ions  along  its  oath;  the  electrons  which  have  a  higher  thermal 
velocity  move  into  this  evacuated  region  before  the  ions  do,  thus  producing  a 
polarized  layer.  Siegel  et  al .  have  proposed  a  mechanism  that  may  give  rise 
to  the  production  of  electromagnetic  radiation  of  about  100  megacycles  as  a 
result  of  oscillations  of  this  polarized  layer. 

The  purpose  of  this  note  is  to  consider  such  processes  and  the  magnitude 
of  the  electromagnetic  radiation  to  be  expected  from  a  satellite  wake.  The 
treatment  to  be  used  here  is  approximate,  but  it  should  be  satisfactory  to 
indicate  the  value  of  this  mechanism  for  detecting  satellites  by  the  radio 
signals  that  their  wakes  generate. 

We  will  assume  for  this  problem  that  a  satellite  moving  7  km/sec  having 
a  diameter  that  will  be  chosen  to  one  quarter  wave  length  (X./4,  about  2  m.) 
passes  through  a  region  with  an  ambient  density  of  10^^  m.  We  will  assume, 
also,  that  the  satellite  can  be  approximated  as  a  flat  disk,  oriented  at  a  right 
angle  to  the  direction  of  motion,  that  simply  sweeps  out  the  ions  along  its 
path,  thereby  producing  an  expanding  (in  length)  cylinder  in  which  the  ions 
have  been  removed.  The  electrons,  because  of  their  much  higher  velocity, 
will  rapidly  move  into  this  region,  thereby,  forming  a  polarized  layer.  This 
layer  persists,  possibly  with  some  oscillation,  for  some  time  thereafter; 
we  will  consider  the  damping  of  these  oscillations  later.  We  assume  that  as 
this  layer  of  electrons  is  formed,  the  electrons  are  collectively  subjected 
to  accelerations  and  that  these  accelerations  of  the  charge  gives  rise  to 
coherent  radiation.  The  steady  polarized  layer,  however,  will  not  give  rise 
to  coherent  radiation  and  therefore  will  be  rather  insignificant. 

At  a  temperature  of  1500°K  the  velocity  v^  of  the  electrons  is  about 
2.61  X  10®  m/sec.  The  electrons  form  a  layer  one  Oebye  length  in  thickness, 
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■  69  ^T/Nq' in  meters  or  about  8.45  x  10  ®  m  for  the  present  case;  thus  their 

acceleration  a  is 

a  =  V  */2\-  ^  4.01  X  10^  m  sec  ® 
e  0 

The  time  required  for  the  electrons  to  move  a  distance  is: 

t  -  2\„/v 
0  e 

This  time  corresponds  to  1/4  cycle  of  the  radio  emission.  The  frequency  f 


and  the  wavelength  \  are  therefore: 

f  =  Vg  /Nq/T'  552  =  3.84  Me 
\  ■  Y8  m 

The  total  charge  participating  in  this  motion  at  any  instant  is  given  by  the 
amount  of  charge  lying  in  a  ring  about  the  satellite  of  radius  thickness 
which  is  the  thickness  of  the  polarized  layer,  and  length/^.  The  radius 
R  has  been  chosen  to  be  \/8  since  it  will  allow  the  opposite  sides  to  emit 
radiation  that  will  be  in  phase.  The  length  of  the  ring  is  approximately 
given  by  the  distance  traveled  by  the  satellite  in  a  time  'K  . 

i  -  v^/v^ 

where  is  the  satellite  velocity.  Thus,  including  an  extra  factor  of  1/2 
to  take  care  of  the  results  which  would  be  obtained  by  an  integration  of  the 
contributions  to  the  amplitude  of  the  electromagnetic  wave  from  the  various 
portions  of  the  ring,  we  get  for  the  total  number  of  electrons  N: 

The  total  power  is  given  by 

_  e^  a^  TTe^  ..  z  2  2 


bn  60  c' 


v^*  v/ 


where  e  is  the  electronic  charge  and  c  is  the  velocity  of  light.  If  we 


substitute  into  this  expression  the  values  for  X.  and  as  they  depend  on 
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T  and  No  we  have 


P  - 


TTe*  1*  Vs 


(69)‘ 


6  Co  c 

in  MKS  units.  This  is  independent  of  both  the  frequency  of  the  radiation 
and  of  the  ambient  density  of  the  plasma.  For  our  case  this  yeiids 
P  -  1 ,26  X  lo"*'"*  watts 

The  amount  of  power  that  falls  on  the  receiving  antenna  is  reduced  by 

the  factor  r*/4d*  where  r  is  the  antenna  radius  and  d  is  the  distance  to 

the  satellite.  Thus,  the  power  falling  on  a  60  foot  antenna,  p  ,  when  the 

d 

satellite  Is  100  miles  away  will  be  reduced  by  a  factor  of  8.1  x  lO**-®. 

Thus,  p  v/ill  be  1.02  x  10  watts. 

Apparently,  the  smallest  amount  of  power  on  the  antenna  which  can  be 
detected  by  the  60  foot  radio  telescope  is  about  1.5  x  lo"®®  M  (watts 
where  Af  is  the  bandwidth.  If  Af  is  taken  to  be  l°/o  of  the  frequency  of 
radiation,  i .e.  56  Kc,  then  the  minimum  detectable  total  power  would  be 
about  5.4  X  lo"*’®  watts.  The  galactic  noise  level  is  about  lo’^®  watts/til® 
at  Me  which  is  about  six  orders  of  magnitude  greater  than  the  power  from 

the  wake.  Thus  it  v.ould  be  possible  to  detect  such  radiations  only  from  the 

satel i i lc  i tself . 

If  we  now  return  to  the  question  of  the  damping  of  these  oscillations, 
we  can  obtain  an  idea  about  the  total  electromagnetic  radiation  that  is 
produced  in  the  entire  wake.  We  merely  need  to  note  that  the  Landau  damping 
formula  shows  that  for  plasma  wave  oscillations  with  a  wave  length  shorter 
then  X.JJ,  the  wave  is  damped  out  rapidly  while  for  those  oscillations  with  a 
wavelength  much  greater  than  the  oscillation  will  be  damped  out  only  very 
slowly.  Since  the  wavelength  for  the  plasma  oscillation,  is  about  ^.jj/10 
for  the  case  considered  here,  the  oscillation  will  be  damped  out  immediately, 
thus  agreeing  with  the  treatment  given.  In  any  case,  however,  we  could  not 
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expect  a  situation  to  arise  where  the  osciilations  would  continue  for  some¬ 
thing  like  10®  X./\  oscillations  as  would  be  necessary  to  make  the  signal 
P 

observable  with  the  60  foot  radio  telescope. 

Tliere  is  also  in  this  approximate  treatment  some  degree  of  error 
produced  by  the  artificial  removal  of  the  satellite  from  the  problem.  In 
the  actual  case,  the  satellite  is  still  very  much  in  the  neighborhood  when 
the  electrons  become  subject  to  accelerations  by  the  electric  field.  It  Is 
possible  that  the  entire  front  surface  of  the  satellite  contributes  to  the 
production  of  radiation.  Nevertheless,  because  of  the  low  velocity  of  the 
satellite  compared  to  that  of  the  electrons,  it  does  not  seem  that  this 
would  significantly  change  the  results. 

Although  the  detection  of  such  radiation  appears  to  be  out  of  reach 
for  the  radio  telescope  at  the  present  time,  it  seems  that  this  radiation 
could  be  easily  detected  from  the  satellite  itself.  The  author  proposes, 
therefore,  that  an  experiment  be  performed  on  an  artificial  earth  satellite 
to  try  to  observe  the  electromagnetic  radiation  over  a  wavelength  of  frequency 
that  may  be  produced  in  the  satellite's  own  wake. 


